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SOME OPERATING- TERMS PERTAINING TO COOLING TONERS 

App r oach : The difference in degrees centigrade between 
the temperature of the cold water leaving the cooling 
tower and the ambient air wet bulb temperature. 

CeU.: The smallest tower subdivision which can function 
as" an independent unit. Each cell he,s its own air and 
water flow system with one or more fans or stanks. 
C oolin g ra nge ; The number of .degrees in centigrade, 
water is cooled in the tower. It is the difference 
between the inlet hot water temperature and the outlet 
water temperature. 

Drift j The loss of water in the form of fine droplets 
being carried away by the exhaust air. Drift is measu- 
red as the percentage of the circulating water. 

Evapo rn,tion~ l oss : The amount of water evaporated from 
the circulating water into the atmosphere. Expressed 
as the percentage of total water flow rate. 

Do g gin g: ¥hen the warm saturated air is discharged out 
of the tower, it comes in contact with the cclder at- 
mospheric air resulting in the formation of fog. 

Heat l oad : The amount cf he.at dissipated in a cooling 
tower in kcal/min from the circulating water. 

Make up : The quantity of water added to the basin to 
replace water lost by evaporation, drift, blow down 
and leakage ( if any) . 

Recircula tion: The entry of a portion of the discharg- 
ed air along with atmospheric air into the air inlet. 
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ABSTRACT 

A cooling tower is an enclosed device for the 
evaporative cooling of water by contact with the air. 
This is achieved partly by an exchange of latent heat 
resulting from the evaporation of some of the circula- 
ting water, and partly by a transfer of sensible heat. 

Cooling tower industry has a very competitive 
market, and hence the refinements in this field have 
been considered of primary importance. The manufactu- 
rer is required to have a set of guaranteed performan- 
ce curves to refer in selecting a cooling tower for a 
particular application under specified conditions. 

These cover the types of tower and packing which are 
carried by an individual firm. These also cover the 
operating conditions such as water flow, cooling range, 
cold water temperature, wet bulb temperature etc., in 
order to design a cooling tower. 

The cooling tower manufacturers in our country do 
not have guaranteed performance curves and have to, 
therefore, either guess the tower size or they depend 
on their foreign collaborators. The present work has 
been undertaken as an attempt towards • solving this 
problem of the cooling tower industry. Performance 
curves have been drawn for both the counter flow and 
cross flow mechanical draft cooling towers for power 
plants, fertilizer and air conditioning plants, de- 
signed to be located in big industrial cities of India. 
Generalised computer programs have been developed, 
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based upon the cooling tower theories already developed. 
Results have been analysed and discussed. 

It is expected that the performance curves deve- 
loped in the present work should be of great help to 
the cooling tower industry in the country and to the 
buyers in selecting and predicting tower performance 
at varying operating conditions. 



CHAPTER I 


INTRODUCTION 

Conservation and reuse of processed water have 
become a necessity world over. Large industrial 
water users are power generation plants, chemical 
plants, steel plants, petroleum refinaries, atomic 
power plants, air-conditioning and refrigeration 
industry etc. An immense quantity of water is used 
by most of the above industries for cooling purposes, 
for example, in condensing the exhaust steam in power 
plants, in liquefying the chemical products in vapour 
state and, in many cases, in preventing overheating 
of the machinery parts which are exposed to high 
temperatures. The standards required of the cooling 
water as regards its temperature and quality, i,e. 
its contents of impurities, may vary considerably 
depending on the purpose of the cooling. water. 

It is necessary that the temperature of cooling 
water should not exceed a certain prescribed value for 
a particular process plant and that its content of 
impurities should not result in the formation of 
deposits in the system or in corroding the metal parts. 
These requirements are dictated by the nature of the 
production processes and by the need for reliable and 
economical operation of the plant concerned. A rise 
in temperature of the cooling water used in a steam 
plant, for example, increases the fuel consumption in 
power production and lowers the plant* s capacity j in 





refrigeration, it effects the coefficient of perfor- 
mance of the plant and in oil refinaries and chemical 
plants, it lowers the yield of the products. Similar 
effects are experienced because of the impurities 
present in the cooling water. Another of the more 
important conditions on the temperature of the cooling 
water is that it should not be heated to a very high 
temperature in the condensing plant i.e. large quanti- 
ties of heat should be transferred at a law temperature 
to achieve the most efficient results. This 
necessitates a very high consumption and continuous 
supply of fresh water. 

Due to the increase in power demand, expansion of 
industries, etc., the use of water has more than 
doubled in the past decade, and its resources, every- 
where, are limited. One has to mainly depend upon 
seas, rivers, lakes, ponds and underground systems as 
sources of water supply. In tropical countries like 
India even these resources are not easily available in 
most of the parts. Underground reserves [l] are 
becoming exhausted because of heavy exploitation over 
a considerable period and a lack of legal control* 

Thus, at many places, rain water is the only source 
left. 


When a sufficiently large supply of water is 
available from the sea, river, lake etc., a continuous 
water supply system is mostly used, in which the water 
taken from the source is used once for the cooling 
purposes and is then discarded. Considerable variations 
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in the water level in a river or lake are found and 
sometimes the water is to he transported over large 
distancos or to a great height. Extensive mineral- 
ization of water or contamination by chemically - 
aggressive impurities, requires an expensive continuous 
purification process. Under such circumstances, when 
the procurement of usable water from its source is so 
expensive, cooling by a continuous fresh water supply 
cannot be adopted. 

Problem is not only to obtain a continuous fresh 
supply of water in a sufficient quantity from the 
source for cooling purposes in industries, but an 
equally important problem is to find a way of discard- 
ing water that has been used for cooling. Temperature 
is a primary factor [2] in the solubility of atmo- 
spheric oxygen in water, and in all chemical and 
biological processes occuring in water. The discharge 
of large quantities of heat to receiving waters may 
affect important processes adversely and result in a 
less desirable quality of water. Of all water uses, 
the propagation of aquatic life is perhaps the most 
affected by temperature. For this reason, establish- 
ing temperature standards to protect fish and their 
food organisms has been of primary importance. 

Reuse of industrial cooling water is a big step 
not only in the water conservation program, but it 
also eliminates the problems of discarding hot water 
resulting in thermal pollution and is very economical 
too. 



4 


One of the methods to cool water is to use air as 
an external heat absorbing medium when hot water is 
being circulated through the tubes. Water acts as an 
intermediate heat - carrier between the plant being 
cooled and the outside air. The heat removal takes 
place by conduction through the tube wall and by 
convection from the outer surface of the tube to the 
moving air. 

The use of air as an external heat — absorbing 
medium has not been adopted very widely. The average 
value of the convective heat transfer coefficient of 
air varies from 2.7142 x 10”^ kcal/sec - m 2 - °K (2 Btu/ 
hr - ft 2 - °P) to 13.571 x 10“^ kcal/sec - m 2 - °K 
(10 Btu/ hr - ft 2 - °F) depending upon the velocity and 
temperature of air [1]. Hence, due to the poor heat 
transfer from the surface being cooled to the air, the 
cooling surface required is many times greater than 
that required for cooling by water. Also, the specific 
heat of air is low and, hence, a high power consumption 
is required for fans to supply large quantities of air. 
Nor it is possible with this device to cool water 
below the ambient temperature of air. Therefore, this 
device cannot be used as an inexpensive and effective 
device for cooling water in industrial units. 

Considerable increase in the rate of heat trans- 
fer between the atmospheric air and the circulating 
water can be achieved by bringing hot water into direct 
contact with moving air. This employs the principle 
of evaporative cooling of the water. In evaporative 
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cooling, a part of the water gets evaporated and the 
water vapour thus produced carries away with it heat 
which is called the latent heat of vaporization. The 
evaporative cooling approach is basically a water 
saving technique. 

The thermodynamic principle of evaporative cooling 
is that water must have heat to change from the liquid 
to the vapour state and when water evaporates, this 
heat is removed from the water remaining in the liquid 
state, dropping its temperature. It takes approxi- 
mately 560 kcal (1000 Btu) to evaporate 1 kg ( 1 lb) 
of water. Therefore, the available cooling effect for 
an evaporative type cooler is about 560 kcal per kg of 
water evaporated. 

In a once - through system using city, well, or 
surface water, each kg of water circulated will pick 
up only 1 kcal for each 1°C of temperature rise, since 
only 1 ’sensible' ' heat gain through the unit being 
cooled is involved. If the temperature requirements 
of the heat source allow for a 10°G temperature risg, 
the heat pickup would be 10 kcal/ kg of water 
circulated. Therefore, in an evaporative cooling 
device 1 kg of water does the same work as 56 kg of 
water in a once-throtigh system. Theoretically, 
evaporative cooling requires only about 2 per cent of 
the water compared to a once -through system. The 
removal by air of latent heat plus sensible heat, 
makes the water cool by evaporation in a cooling 
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device, as shown in figure (1.1), 


f f 


pump 



cold water 

Fig. (1.1), Schematic diagram of circulating cooling 
water-supply system. 


The air is chosen as a tiLedium for cooling water 
because it is so freely available in -unlimited 
quantities* It not* only absorbs the water waste heat 
load but also carries it away and dissipates it into 
the atmosphere without altSKiz® the atmospheric 
conditions. When evaporative cooling is used, the 
cost of water collection from the source and its 
purification is incurred only once and the waste water 
problem is almost fully eliminated. 

Because of these advantages, evaporative cooling 
of circulating water has predominated in recirculation 
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water •• supply systems. 

Figure (l.Sj shows a portion of a psychrometric 
chart and indicates the psychrometric analysis of the 
air path through the evaporative cooler. The true 
path is approximated by the curved dotted line from 
point A (entering air conditions) to the point C 
(leaving 'air conditions) . This actual path will vary 
either above or below that shown, depending on the 
unit design and the ratio of air to water flow. 



Fig. (1.2), Cooling effect of air passing through 
evaporative cooler. 

For purposes of explanation, the air path is 
broken down into two vectors, lines AB and BC. Air 
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enters the evaporative unia in an unsaturated ambient 
condition (point A) . Before reaching the heat transfer 
surface, it is saturated without change of total heat 
content as it travels to the point B. Passing across 
the heat transfer surface, the air absorbs heat from 
falling water. As some of the rater evaporates * the 
heat content of the air is increased as the temperature 
of the falling water is reduced. Since the air is 
continually being washed by the rater, the process 
follows the saturation line to the final leaving air 
temperature (point C) . Travel from point A to point B 
is adiabatic there is no cooling of the water. There 
is only a conversion of air sensible heat to latent 
heat as the air temperature drops 'oo that of the wet 
bulb temperature. The effective heat removal takes 
"'lace between the points B and G along the saturation 
line. Thus, the wet bulb condition of the entering 
air at the point A is the only factor affecting the 
cooler performance. 

The principle of evaporative cooling was well 
known since ancient times in India and other Asian 
countries. Only recently, however, it has been put to 
practical industrial application in these countries. 

In U.S.S.R., U.S.A., and European countries where 
industrial development started much earlier than in 
Asian and African countries, evaporative cooling 
gained a great deal of importance and has been long 
analysed and studied scientifically. 

Among the earliest theoretical and experimental 


f 
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studies made in England on this subject, I,V. Robinson's 
work in 1907 [3] on the theory of cooling towers appears 
to be the first in the field. Among the earliest 

> 

studies made in the U.S.S.R. , were those of A. IT. Arefev 
in 1925 at the E, Dzerzhinskii All - Union Power 
Engineering Institute on cooling towers [4] and those 
of N.M. Bemalskii in 1929 on cooling ponds [4] . These 
studies formed the basic of the development of the 
theory and methods for the thermal design of coolers 
employing evaporative cooling of water. It was only 
in 1926, that a successful theoretical analysis of 
heat and mass transfer in cooling towers was made by 
Merkel in Germany [3]. The Carrier charts for cal- 
culation of moist air relationships were developed in 
U.S.A. in 1911 and a similar chart by Mollier in 
Europe in 1923. Extensive efforts to explore the 
field of industrial cooling plants were not made until 
after the incidence of second world war. Design of 
various types of cooling devices and their use for 
industries, homes and buildings become familiar all 
over the world around the 1930 * s . 

The main types of water cooling devices are 
classified as follows: 

1) Ponds 

(a) Cooling 

(b) Spray 

2) Closed Coolers (natural draft) 

(a) Spray - filled 

(b) Splash - filled and 
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3) Tower Coolers 

(a) Natural draft chimney tower 

(b) Mechanical draft tower 

(i) Forced . draft 

(ii) Induced draft (counter- flow, 
cross-flow) 

1(a) Cooling Ponds : It is the cheapest and the simplest 
method of cooling water. It consists of a comparatively 
large pond in which cooling takes place by air contact 
at the surface. In order to achieve the best results, 
the total area of the water surface should be as large 
as possible. This is achieved by having the inlet and 
outlet points placed as far apart as possible and by 
raising the water level in the pond. 


The equation for estimating the size of an evapo- 
ration pond needed for a specific incoming flow rate 
is [5]s 


D 


( 


325.600 L 


7.48 x 43,560 x A 


E 

n 

12 




( 1 * 1 ) 


where 


D = epth of pond, ft. 

L = Flow of processed water into the 
pond, gpm. 

A = ;rea of pond, acres. 

E n = Net evaporation (evaporation - rain 
fall) , inches/ year. 

T = Time, years. 


conversion factors: 

525,600 = mins per year 
7.48 = gal/ ft 3 
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43,560 = ft 2 / acre 

12 =' in/ ft 

A cooling pond has the following advantages % 

(i) Cooling pond may be constructed at a very 
low cost by pushing up a certain height of 
the earth. 

(ii) It may operate for a long period without 
requiring any make-up water and with low 
maintenance cost. 

However, the use of the cooling ponds is quite 
limited because of the following disadvantages i 

(i) The heat transfer rate from a cooling pond 
is very low. The heat dissipated from a 
still pond averages 17 kcal per hour per 
square meter of water surface of a still 
pond, per degree centigrade temperature 
difference between the contacting air and 
water surface. 

(ii) large areas needed for the ponds create 
serious problems particularly in big indus- 
trial cities. 

1(b) Spray Pond : A basic modification of the cooling 
pond concept involves the use of spray ponds. The 
spray pond cooling system pumps heated water through 
spray nozzles, which divide the water into small drop- 
lets, thus increasing the effective area for evapora- 
tive heat transfer to the atmosphere. Cooling occurs 
in the spray pond as the water is propelled upward and 
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falls to the surface of the pond* The pond acts largely 
as a collecting basin. 

The amount of cooling to be gained in direct air- 
water contact of the spray system ultimately depends 
on the air-temperature, humidity and wind conditions. 
Higher wind speeds effect more efficient heat transfer 
to the atmosphere. A lou/red fence may be provided to 
reduce the loss of water with the outgoing air. 

Although spray ponds are more compact in design 
and better in their performance than cooling ponds, they 
still require larger areas and have limited performance 
capacity because the time of contact of sprayed water 
and air is too small. Their use, therefore, has nearly 
completely stopped since the cooling towers came into 
existence in the 1920’s. 

2(a) Spray-filled Cooler (Fig 1.3): 


hot 

water 

inlet 


air 


cold < 

water 

outlet 



air 


Fig. (1.3), Spray-filled cooler 
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This is essentially a narrow small-sized spray tank 
with elevated nozzles and a high louvre fence. The air 
movement through the louvre d. side is fully dependent 
upon wind conditions. Since wind normally blows in a 
horizontal direction, flow of air is crosswise to the 
flow of water. 

The spray -nozzles used are of the tangential type 
which are less liable to become clogged up. The operat- 
ing pressure before the nozzles is from 3 to 5 meters 
of water, or sometimes, between 12-14 meters of water. 
Nozzles pointing downwards are used which results in 
more uniform distribution of the water and reduces drift. 

Spray-filled coolers are suited for refrigeration 
and air-conditioning system applications. Table (1.1) 
gives the typical dimensions of spray-filled coolers, 
with nozzles pointing downwards [4]. 

These towers can be installed in the open where 
there is no obstruction to prevailing winds. It can 
be installed on the roof of a building or on a special 
platform raised for the purpose, if there is a short- 
age of space. 

A spray-filled tower has the following disadvant- 
ages; 

(i) The approach to the wet bulb temperature 
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Table (1.1) Typical dimensions of spray coolers, 


1.1 
2.3 
4.5 
10.0 
20.0 
30.0 
50.0 
50 .0 
50.0 
100.0 
100.0 
100.0 
150.0 
200.0 

250 .0 

300.0 

340.0 


Area 

2 


0.85 

1.50 

1.90 

3.40 

4.60 

8.60 

20.60 

21.40 
21.00 
36.80 
38.80 
46.30 
61.50 

76.40 
91.70 



width 


920 
1220 
1220 
1830 
2140 
2140 
2140 
2750 
4020 
2750 
4020 
4020 
4020 
4020 
40 20 
4020 
4020 


Dimensions, mm 


length ? Height 


920 

1830 

1220 

1830 

1550 

2750 

1830 

2750 

2140 

3660 

4020 

3050 

9660 

2750 

7780 

2750 

5240 | 

2750 

13400 

3660 

9660 

3660 

11540 

3660 

15300 

3660 

19000 

3660 

22800 

! 3660 

26600 

3660 

30500 

3660 





will always be equal to or greater than the 
cooling range, except when relatively high 
hot water temperatures are encountered. 

(ii) There are windage losses. 

2(b) Splash - filled Coolers ; The main difference 
between the splash - filled coolers and spray - filled 
coolers, is that in the splash - filled coolers, fill- 
ing is used to increase water break-up and to provide 
additional water surface to air flow. The cooling 
efficiency is greater in splash - filled coolers. It’s 
increased cost and maintenance makes it obsolete. 


3(a) Natural - draft Towers : Natural- -draft towers were 
the first large cooling apparatus built. It consists 
of an empty shell made of steel - reinforced concrete 
structure mostly in hyperbolic designs. These are 
built very high up to 120 meters with the base 
diameter 80 meters. At the lower end is the packing 
through which water trickles, drops or flows in a 
predetermined manner so as to give up a portion of 
its heat to the air stream flowing past. 

In the hyperbolic-type tower, shown in figure 
(1.4) the packing occupies the whole of the space 
between the hot water distribution system and the air 
inlet position. There is no air distribution system, 
and' the air enters the packing at the periphery of 
the bottom of the tower. 
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Fig. (1.4), Hyperbolic Tower. 

The hyperbolic shape of the tower does not have 


cold water 
out 
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any thermodynamic significance. It could be of a cylin- 
drical shape as well. 

From an aerodynamic point of view, hyperboloids 
or cones are better than cylinders [6]. Hyperbolic 
tower is still better than the conical as it directs 
the entering air more smoothly towards the centre, and 
the upper rim tends to produce a stronger upward draft 
than the inclined straight rim of the conical. The 
hyperbolic tower, since it is a doubly curved shell 
has a great strength and shearing stresses are elimi- 
nated. 


Air flow in the hyperbolic tower is produced due 
to so many factors. The difference in the specific 
weights of the cold dry air outside and the air inside 
the chimney # and the temperature and humidity raised 
because of its contact with warm water are responsible 
for the air-flow movement in the hyperbolic tower. 

The most significant contribution is due to the draft 
created by the height of the chimney. The greater the 
height, the greater is the draft achieved. The ave- 
rage velocity of air above the packing is of the order 
of 1.5 - 2 meters/second. 

The main advantages of natural draft cooling 
towers are the followings 

(i) No fans are required to blow air in the 
tower. This not only eliminates the 
capital investment for the mechanical 
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equipment and the related electric control, 
in comparison with mechanical draft towers, 
but it also greatly reduces the expenditure 
towards the operation and maintenance costs. 

(ii) They can practically never break down. 

(iii) They can cope with tremendous water loads. 

(iv) G-round fogging and recirculation of warm 
air are practically avoided because of the 
large height of the tower. 

(v) Loss of water due to drift is negligible. 

The principal disadvantages are: 

(i) The great height necessary to produce the 
draft. 

(ii) Inlet hot waiter temperature must be kept 

hotter than the air dry bulb temperature. 

(iii) Exact control of outlet cold water temp- 
erature is difficult to achieve. 

3(b) Mechanical - draft Towers : Mechanical -draft cool- 
ing towers have a relatively simple construction, 
although the heat transfer processes which take place 
within them are extremely complex. Except for ambient 
air wet bulb temperature, they are virtually indepen- 
dent of atmospheric conditions. They can cool water 
to temperatures below ambient with low capital and 
running costs in a comparatively small space. 

There arc various types of mechanical - draft 
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cooling towers. Figure (1.5) illustrates the ma in 
components common to most towers. 

(i) Casing 2 This is a structure which encloses 
the heat transfer process and provides a 
support for the other main items. 

(ii) Fan : In mechanical -draft towers, a fan is 
fitted to move the required amount of air 
through the water to be cooled. 

(iii) Drift Eliminators 2 These are placed at or 
near the air outlet, and prevent droplets 

of water from being carried out of the tower 
by the air stream. 

(iv) Water Distribution System 2 For maximum 
effect the water entering the tower must 
be spread evenly over the top of the pack- 
ing. Nozzles are used to atomise the water, 
for a spray distribution system and trough 
or weir, where the water spreads by gravity. 

(v) Packing 2 This provides a large water 
surface area to assist heat transfer. This 
may be either splash packing or film pack- 
ing. The depth of the packing may be as 
much as 7 to 8 meters. Th ; two types of film 
packing, now in general use are 2 

(a) grid packing 

(b) plate packing. This was developed to 
increase still further the heat transfer 
surface per unit volume. It consists of 
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Figure (1.5), Main components of a mechanical draft cooling 
tower. 





closely sp-'cocl vertical plates, each 
plate "being corrugated in some manner. 


Most splash and grid packings are of timber 
(specially - treated redwood) , although 
they are also made with T-shaped plastic 
sections . 

Film plate packings have been made in 
galvanized mild steel, anodised aluminium, 
stainless steel, cement, asbestos and 
various plastics [7] . 

The mechanical -dr aft cooling tower uses fan power 
to push the outside air into the tower if it is a 
forced draft or to exhaust warm air outside, if it is 
an induced draft. 

Forced draft : In forced draft towers, the fans and 
the electric motors are easily accessible for mainten- 
ance. Vibration is kept down because the mechanical 
equipment is near the ground and on a solid foundation. 
Some of the disadvantages are, because of the low 
height of these towers. The chances of the hot, humid 
exhaust vapours coming out of the top of the tower, to 
be recirculated are nearly cent percent. The size of 
the fan is also limited to about 4 meters or less to 
suit the tower design, which implies that installation 
of a larger number of towers can meet the demand. 

These disadvantages have been overcome in induced 
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draft towers and, hence, induced draft towers have 
gained preference and special attention by designers, 
sellers and users in the recent past. 

Induced draft (figure 1 , 5 ) • This has the fan located 
"t the top of the tower to push out the warm air of the 
tower. 


This has been further sub-divided ass 

(i) Counter-flow. 

(ii) Cross-flow. 

In counter-flow type, air is drawn up vertically 
against the water falling down in the opposite direc- 
tion through- the packing. Maximum performance is 
achieved in this tower because the coldest water comes 
into contact with driest air and the warmest water 
contacts the most humid air, thereby, having the 
maximum enthalpy gradient at the top at all the sec- 
tions , 


In the cross-flow type, air enters the packing 
horizontally throughout its height through the lomras 
while water is falling down across the packing. 

1 Both the types of towers have advantages and 
shortcomings . 

from a thermodynamic standpoint, the counter- 
flow type is more efficient than the cross-flow type [8] 
because its enthalpy potential difference is higher. 
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Lower static pressure is encountered in the cross-flow 
tower because not all the air passes through the 
entire fill section reducing operating and maintenance 
costs. The cost of both the types is about the same. 

The cross-flow arrangement has an advantage over 
the counter-flow system that the drift losses are less* 
resulting in lower horse power requirements. 

The cross-flow design has a larger plenum area [9], 
and is more susceptible to biological attack. Its 
open-pan distribution system allows the development of 
slime and algae, which must be controlled with chemi- 
cals, but cleaning the nozzles in a cross-flow tower 
is simple. The cross-flow tower can be easily altered 
to increase water flow. 

The counter-flow tower has a closed-pipe distri- 
bution system to control algae, but is harder to clean. 

The co u nter— flow tower performs better than the 
cross-flow type with fans off because of the natural 
draft produced by fewer air intakes. The natural draft 
is particularly important in the winter or in the 
evening when more desirable psy chrome trie conditions 
exist for the tower. 

In India, the first cooling tower was manufactur- 
ed for the Ahmedabad Electricity Company by Gammon 
India Ltd., Bombay around thirties. Prefabricated 
cooling towers for the plants like chemical industry 
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and for central air conditioning plants continued to 
be imported in the country till recently. It was only 
in the last decade, that the following concerns [1] 
started manufacturing over 400 big and small complete 
range of cooling towers per annum in the country with 
different foreign collaborations. 

fable (1.2) shows the main cooling tower manufactur- 
ers in India and their collaborators. 


Table (l.2^ Main cooling tower manufacturers in India 
and their collaborators. 


Si. No. 

Manufacturer 

Collaborator 

1 

Air Conditioning Corp- 

Heenan Proude, U.K. 


oration Ltd. , 

i 


17 Taratalla Road, 



Calcutta - 33 . 


2 

Gammon India Ltd., 

L.G. Mouchell and 


Prabhadevi, Cadell Road,) Partners, U.K.(till 

** \ 

Bombay - 28. 

: recently) 

3 

Larsen and Toubro Ltd., 

Film Cooling Towers 

1 

Dougall Road, 

Brentford, 


Ballard Estate, 

Middlesex, 


Bombay - 1. 

U.K. 

4 

Paharpur Cooling Towers 

The Marley Company, 


Pvt. Ltd., 

Kansas City, 


1-B Judges Court Road, 

U.S.A. (till recent 

» 

Calcutta - 27. 

..... . ■ 

past) 
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Table (1.3) shows some of the important cooling 
tower installations in the country. 


Table (1.3) Important cooling tower installations in India 


Sl.Ho. 

, Cooling tower 

installations 

Capacity of 
Cooling, u 3 /hr 

• Service 

! 

j 

1 

Thermal Power Plant, 
Panki, Kanpur. 

16,000 

Power Plant 
i 

1 

2 

Durgapur Projects Ltd., 

42,000 (3 units) 

> > ; 


Durgapur, W.B. 

' 

1 

3 

i Talcher Thermal 

37,140 

" ! 

| 


; Scheme, Talcher, 
Orissa. 


1 

i 

i 

4 

Ahmedabad Electricity 

26,450 (nat. draft, 

t ! 


Supply Co., Sabarmati. 

15 units) 


5 

Andhra Pradesh 

21,400 (nat. draft. 

t l 

j 

1 Electricity Baard, 

3 units) 



' Kathagudam. 

i 

! 


6 

: Renusagar Power Supply, 
Renusagar. 

14,000 j 

< 

i 

i 

| 

f ! 

7 ■ 

Hindustan Steel Ltd., 
Rourkela, Orissa. 

50,000 

Stool Plant 

8 

Durgapur Steel Plant, 
Durgapur, W.B. j 

19,550 (nat. draft) 

f 

I 

1 1 

9 

Fertilizer Corpn.of j 

India Ltd.,G-orakhpur . 1 

|2l,250 (4 units) I 

i i 

1 - i 

i S 

Fertilizer 

10 

i 

Neyvelli Lignate Corpn|80,000 | 

Heyvelli, Madras. j 

fertilizer 

and Power 


PI ant 
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■ 11 

Gujarat State 
Fertilizer Co., 

Bar o da. 

i 13,000 

1 

\ 

t 

Fertilizer 

\ 

j 

12 

Indian Explosive Ltd., 
Panki, Kanpur. 

13,800 (3 units) 

1 1 

13 

Synthetics and Chem- 
icals ltd. ,Bareille, 

U.P. 

8,500 

Synthetic 

rubber 

14 

National Organic 
Chemicals Corpn. ltd., 
Bombay. 

14,000 

i 

! 

I 

! 

Petrochemic- 

als 

15 

Indian Oil Corpn. 
ltd., Kovali, 

; Gujarat. 

13,316 

Oil-refiner 

16 

Worli Dairy Scheme, 
Bombay. 

1,363 

Dairy 

17 

Hindustan Photofilm 

Mfg. Co. ltd., 

Ootacamund. j 

L 

3,480 

Photofilm 


Many improvements have been made by the individ- 
ual cooling tower companies on the basic theories used 
for analysis. Because the industry has a very compe- 
titive market, these refinements have been considered 
proprietory and closely guarded. Several papers have 
been written by the companies describing the theory 
but few publications are available concerning the 
practical applications , 

, , fhe manufacturer is required to have a set of 
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guaranteed performance curves covering operating con- 
ditions such as water flow, cooling range, cold water 
temperature, and wet bulb temperature etc., in order 
to design a cooling tower. The cooling tower manufac- 
turers in our country do not have guaranteed perfor- 
mance curves which are related to Indian conditions. 
These firms mainly depend upon their collaborators. 

In the present work, known and proven theories 
are applied and computer programming has been develop- 
ed for computing guaranteed performance curves for 
both counter flow and cross flow mechanical draft 
cooling towers for power plants, fertilizer and air 
conditioning plants for Indian conditions. Various 
major cities are considered in studying the performance 
curves of cooling towers with varying wet bulb temp- 
eratures since wot bulb condition of the entering air 
is one of the most important factors affecting the 
cooler performance. 

The curves developed in the present work should 
be of great help to the cooling tower industry in the 
country and to the buyers in selecting and predicting 
tower performance at varying operating conditions. 



CHAPTER II 


FUNDAMENTALS OP COOLING TOUER ANALYSIS 

The generally accepted concept of cooling tower 
performance was developed by Merkel [3] in 1925. His 
analysis and equations include the sensible and latent 
heat transfer into an overall heat and mass transfer 
process based on enthalpy difference as the basic 
driving force. 

Consider a counter-flow tower of 1 m 2 ground 
area through which G kg/hr of air is flowing upward 
and L k^hr of water is flowing downward. The counter- 
flow tower can be resolved into a one-dimensional 
problem [10] with the assumption that the flow pattern 
is vertical with the water falling downward through 
the tower and the air being forced upward. 

Each water particle is surrounded by a film of 
saturated air at the bulk water temperature as shown 
in figure (2.1). The air is being heated and satu- 
rated as it passes through the tower. The heat is 
transferred from the interface to the main air mass by 

(i) a transfer of sensible heat and 

(ii) by the latent heat equivalent of the mass 
transfer resulting from the evaporation of 
a portion of the bulk water. The two 
processes are combined, into a single 
equation 

L dT = K a dv (h* '-h) = G dh (2.1) 
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where, 

L = Vater flow rate, kg/lir. m 2 of tower cross- 
section. 

dT =- Temperature differential, °C . 

K = Overall mass transfer coefficient, 

kg/ hr (m 2 of contact area) (kg water/ 
kg dry air) . 

a = Interfacial contact surface, m a /m 3 of 
tower volume. 

dv = Differential volume, m 3 . 
h ,( = Enthalpy of saturated air at water 
temperature, kcal/kg dry air. 
h = Enthalpy of main air stream, kcal/kg dry 

air. 

0 = Air rate, kg/hr. m of tower cross-section. 
V = Tower volume, m 3 /m 2 of tower cross, 
section. 



Pig. (2.1), Heat and Mass transfer between water, water 
vapour-film and air. 
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In equation (2.1), Merkel makes the following 
assumptions: 


(a) The lewis relationship = 1.0 (2.2) 

P 

where, 

a = Coefficient of heat transfer by convection, 
kcal/hr.m 2 °C. 

C = Specific heat of air, kcal/kg °C. 

ir 

which holds true for the cooling tower (in the case of 
water evaporating into air, conductivity is approxi- 
mately equal to diffusivity and hence this ratio is 
unity [ll] ) . 

(b) The air enthalpy can be expressed as 



T + h fg>s W 


(2.3) 


where, 

hfg s = latent heat of vaporisation dbr water, kcal/kg. 
¥ = Humidity ratio of moist air, kg water vapour/ 


kg dry air. 

The equation (2.3) is an approximate one, which is 
derived by neglecting the superheat in the vapour and 
the heat of the liquid corresponding to the vapour 


content. 


Int egrat ing equat ion (2*1) 

r ' 1 

I dT gaV 

h"-h" 1 

T 2 { 
h2 



dh _ KaV 

h* '-h ~ G 


we get 


(2.5) 

(2.5) 
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In arriving at equations (2,5) and. (2,5), the resis- 
tance to mass transfer from bulk water to interface, 
the temperature differential between the bulk water 
and interface and the effect of evaporation have been 
ignored. The left hand sides of the equations contain 
only the thermodynamic conditions for the cooling 
process. It is determined wholly by the initial and 
end conditions of the air flowing through the tower. 
The right hand sides of the equations are independent 
of the thermodynamic conditions in the tower and are 
determined by the characteristic of the toxfer design 
KaV and the water and air flow rates L and G-. 

Equations (2.5) and (2.5)’ are the basic equa- 
tions for calculating cooling tower performance. 

The performance of a cooling tower quite patently 
depends upon a number of variables. The cooling tower 
industry today is called upon to design towers requir- 
ing large ranges and close approaches to the wet bulb. 

In order to construct the performance curves it 
is of course, necessary to choose from among the many 
variables those which are to be used as co-ordinates, 
those to be used as parameters, and those to be held 
constant. This involves the following variables: 

(i) approach to the wet bulb 

(ii) range 

(iii) l/G- ratio 

(iv) wet bulb temperature of the entering air^ 
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The tower characteristic is the required 

factor and, therefore, it must he chosen as the 
ordinate, L/ G- is selected as abscissa, since it is 
one of the most important design characteristics. 

In cooling tower practice, G is usually chosen 
from power requirement considerations. If L/ G is 
known, then L can be obtained and therefore the nece- 
ssary tower ground area for a given tower filling and 
a given tower capacity can be calculated. 

Generally a cooling tower is selected initially 
by reference to sets of performance curves which each 
individual manufacturer has prepared for his own use, 
and which cover the types of tower and packings which 
are carried by that firm. Usually such diagrams 
consist of charts showing the approach to wet bulb 
temperature for different atmospheric conditions and 
varies 'as functions of the tower characteristic^— 

and || . 

Thus for any wet bulb temperature, range and 
approach, and for any chosen ratio, these curves 

give 

£ dT _ 

T h ' 1 -h 

T? J 

Since KaV is a function of L, as well as of G, ' 
for each tower design, curves of as functions of 
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-g should be constructed. 

When analyzing the cooling tower performance, 
the application of the equations become extremely 
complex since, with the cooling tower, if one factor 
of the equation is changed, it starts a chain reaction, 
automatically changing other factors involved. 
Consequently, it becomes difficult to segregate and 
evaluate the effect of each individual factor. 

The performance chart allows to select a cooling 
tower and predict its performance at any other condi- 
tions. It helps in evaluating the performance of an 
installed cooling tower. With one test taken at any 
operating conditions, it is possible to: 

(a) determine if the tower is delivering its 
rated capacity , 

(b) predict the tower's performance at any other 
operating conditions. 

* 

The performance of a cooling tower depends upon 
a number of variables. These are the cooling range, 
the approach to the wet bulb, the entering wet-bulb 
temperature, and the total circulating water flow rate. 
Performance is also affected by the tower and packing - 
design and by the water-to— air ratio required to meet 
the specific design point. 

Packing : The performance of a cooling tower is largely 
dependent upon the packing design. The function of 
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the fill is to increase [12] the contact between air 
and water by offering new exposed surfaces to the air 
as the water flows through the tower. Another func- 
tion is to maintain proper distribution of both air 
end water. 

The performance of a cooling tower is improved 
by increasing the amount of filling, height, area and/ 
or air quantity. Increasing the tower height increases 
the length of time the air is in contact with the 
water, without seriously affecting the fan power 
r equired, but increases the pumping power. Increasing 
the tower area while maintaining the constant fan 
power increases the air quantity somewhat and increases 
the time that this air is in contact with the water 
because of lower velocity. The surface area of water 
in. contact with the air is increased in both cases. 

Considerable work has been reported in the 
literature about the ‘use of various types of packings. 
Lichtenstein [13], London et al [14] and many other 
investigators have reported the performance character- 
istics of wood grids as cooling tower packings. Lowe 
and Christie [15] investigated the packing made with 
both flat and corrugated asbestos cement sheets. 

L aray ankh e dkar et al [16] investigated the packing 
consisted of alternately arranged flat and corrugated 
aluminium strips. This packing gave higher values of 
Ka under approximately same values of L and G com- 
pared to the packings used by other investigators 
[13, 15]. Bulanina et al [17] mainly concentrated on 
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fillings for the high capacity cooling towers rated at 
65,000 - 100,000 m 3 /hr. Aerodynamic and thermal 
oondirtions An cooling towers having a wetted area up 
to 12,000 m 2 were investigated. 3 types of fittings 
were considered. They were flat asbestos-cement panels, 
corrugated asbestos-cement panels, and wooden slots. 

The most suitable arrangement of the packing can 
only be found by experiment. The object is to achieve 
a sufficiently large cooling surface whilst minimizing 
the resistance offered by the packing to the passage 
of air, a loss in the performance as a result of a 
slightly smaller cooling surface con be compensated by 
increasing the air flow and air velocity for the same 
height of stack. Increasing the air velocity through 
the tower decreases the time, the air is in contact 
with the water, but, since a greater quantity is 
passing through, the average differential between the 
water temperature and the wet bulb temperature of the 
air is increased, and this increases the rate of heat 
transfer. Increased air quantities are obtained only 
at the expense of increased fan power, which increases 
approximately as the cube of the fan speed [18] . 

When it is very difficult to determine accurately 
the free surface of the liquid, e.g., on breaking up 
the flow of circulating water into droplets, the volu- 
metric of heat and mass transfer coefficients are 
used, i.e., coefficients that are based, not- on the 
unit surface of the water, but on the unit active vol- 
ume of the cooler [4]. * 
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L.D. Berman [4] found that the volumetric mass 
transfer coefficient Ka for a tower packing is approx- 
imately proportional to the mass velocity of air to 
the power 0.55 - 0.60 and to the superficial water 
flow rate to the power 0.3 - 0.4. These experiments 
are also confirmed by Lichtenstein. 

For a counterflow tower with a packing of rectang- 

* 

)ular slats Lichtenstein [13] obtained an empirical 
relation: 

Ka = A G m L" = B G m l n (2.7) 

The constants in this equation are: 

A = 635, B = 1050 
m = 0.53 and n =0.39. ' 

The size and arrangement of the slat packing are 
given in figure (2.2) . 



Fig. (2.2), Size and arrangement of rectangular slats 
by Lichtenstein. 
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The volumetric mass transfer coefficient for the 
above tower packing as a function of water flow rate 
is given in Table (2.1) . 


Table (2.1) , volumetric mass transfer coefficient for 
the packing in figure (2.2) as a function of water 
flow rate. 


jwater 
flow rate 
l(kg/hr.m 2 ) 



Ea for 

G=5000( kg/hr. m 2 ) 

Ea for 

G=6Q00 (kg/ hr . m 2 ) 

- ■' ■» 

Ea for 

G=8000 (kg/ hum 2 

4,000 

1,050 

1,450 

1,750 

6,000 

1,200 

1,700 

2,050 

8,000 

1,400 

1,950 

- 2,500 

10,000 

1,575 

2,150 

2,530 

12,000 

~ » _ — 

1,700 

! 

2,500 

2,700 


Water Load : For any industrial plant, the method of 
heat dissipation dictates the ^mount of water flow 
through the intake. The amount of water passing 
through the cooling tower, L, is determined by the 
amount of heat to be transferred and the allowable 
increase in the cooling water temperatures 


1 



( 2 . 6 ) 


where , 

Q = Heat removed per unit time, kcal/hr. 

T = Allowable temperature increase of cooling 
water, °C^,. 
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C^ T = Heat capacity of cooling water f kcal/kg.°C. 

For a cooling tower, if the mass air rate remains 
constant, an increase in water rate results in a higher 
wet bulb temperature of the outgoing air, The temp- 
erature of the air leaving the tower can approach, but 
not exceed, the temperature of the water entering the 
tower. 


Investigators have placed an upper and lower limit 
on water loading, some even stating that a loading 
between 5.0 and 12,0 m 3 /hr.m 2 gives optimum [19] per- 
formance, The reason is that a light loading results 
in poor water distribution and unequal wetting of the 
filling while a high loading causes what is referred 
to as flooding. 

Capacity of an existing tower cannot be increased 
by pouring more water over it. Additional water still 
contacts the same amount of air, with the net effect 
of raising the cold water temperature. It is not ad- 
visable to reduce water flow when less output is 
needed. Output should be decreased by shutting down 
cells, reducing fan speed, or turning off fans. 

Evaporative water losses may vary widely, depend- 
ing upon the type of cooling facilities and the exist- 
ing ambient conditions. 

The general water balance assessment of a cooling 
system is expressed as s 
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M = E + D + B (2.8) 

where , 

M = Make up water flow rate, 

E = evaporation rate, 

D = drift and windage loss rate, and 
3 = blowdown flow rate. 

Evaporation loss averages approximately 1 per cent 
of the circulating water for every 5°C cooling range 
[20]. The drift loss on a modern mechanical draft 
tower is considered to be something less than 0.1 per 
cent [2l] of the circulating water flow rate. 

Cooling Range, Approach and WBT ; These affect the tower 
performance in terms of the tower size and cost. The 
broader the range, the more expensive the tower is. 
Longer air-water contact, which brings about greater 
temperature .drops, depends on a larger fill area and 
fan. Thus a broad range will also result in an in- 
creased cost of operation. 

The smaller the approach, the higher the cost of 
the tower. Trying to lower water temperature by 1 or 
2 extra degrees requires dramatic increases in tower 
size and operating costs. 

The design wet bulb temperature should be select- 
ed on economical conditions [22] and will not necess- 
arily be the highest WBT registered in the area. 

As the periods when the actual wet bulb temp- 
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erature exceeds the design WBT are for short duration 
during summer mid-day hours, it is more economical to 
operate at slightly higher water temperature from the 
tower during these hours than to install a tower that 
would he oversize for all hut these short periods. 

It is characteristic for natural draft towers to 
have relatively better output at lower wet bulb and for 
mechanical drafi: towers to have better performance at 
higher wet bulbs. 

Metho ds To Eva l uat e l ower Performance ; 

C PUNIER FLOW TQWEil PERFORMANCE ; The method of correla- 
tion developed by Lichtenstein [13] can be used for 
designing, evaluating and predicting counterflow tower 
performance. Experience and research conducted by the 
cooling tower industry and other organisations have 
shown that this method approximates the performance of 
counterflow cooling towers. The method has gained 
wide acceptance by the industry. 

*1 

r 

| dT 

A graphic representation of l is shown 

in figure (2.3) . 1 pJ 

The water temperature 1 is selected as abscissa 
an d the enthalpy as the ordinate. Saturation line h’ ’ 
gives the enthalpies of saturated air at water temp- 
eratures. In figure (2.3) arc<£ represented by ABCD 
determines the tower characteristic necessary to cool 
water from T-^ to T£ with inlet air enthalpy of h-^ and 

a given water to air flow rate ratio . L/G. The air 
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enthalpy varies as a straight line with the water 
temperature since L dT = G dh. 



Water temperature, T °G 
Fig, (2.3) Graphic representation of 



T, 


ar 

h"-h 


It is possible to introduce a * 1 log-moon potential* * 
as an approximate method of solution in accordance 
with the analogy of the heat exchanger. The log-mean 
enthalpy potential would be mathematically correct if 
the enthalpy potential h* *-h were a straight function 
of T. This is true only if the saturation line is 
straight, and so it follows that the log-mean potent- 
ial can give good results only for very small ranges 
over which the saturation line could, be considered as 
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approximately straight. Test data taken under one set 
of conditions cannot he used accurately to predict 
performance under other conditions if this method is 
employed, since the error changes -with the curvature 
of the saturation line. Por practical cooling tower 
work, the use of a log-mean potential is, therefore, 
not adequate. 


London et al [14] expressed the tower performance 
in terms of its effectiveness as an energy exchanger. 


1 

e h - <3F (h al -h 2 ) 


(2.9) 


where , 

= Unit heat capacity of liquid water, 
kcal/kg°C. 

h a = Enthalpy of the air-water vapour 

mixture at the equilibrium WBT, kcal/kg. 
h = Enthalpy of air stream, kcal/kg. 

Suffix, 1 = Water entrance, air exit. 

Suffix, 2 = Water exit, air entrance. 

If, as normally is the case, the cooling range is 
such that 

S 2> Vb2, L W T 2> ^ G (hi“h 2 ) 

- cK^-hg) 


= h l~ l12 (2.10) 

il al" il 2 

i.e., represents the ratio of actual tower energy 
exchange between phases to the energy exchange which 
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"would result provided the discharged air was saturated 
at the temperature of the entering water. 

may .be represented empirically by an exponen- 
tial equation: 

-n K a V 

e h = 1-0 e & (2.11) 

where, 

C = Dimensionless Coefficient, 
r] = Traction of the packing area aV covered 
by the flowing water. 

In this method, C and t] in equation (2.11) had 
to be found experimentally. The effectiveness was 
found to vary from 0.3 at low water rates and high air 
rates, to 0,8 at high water rates and low air rates. 

Baker et al [19] presented the Unit-Volume Co- 
efficient as a method of cooling tower performance. 

In this method the integration is accomplished by 
increments of equal volumes. The Unit-Volume Coeffi- 

3 

cient can be defined as kcal transferred per m of 
2 

tower per m of plan area per kcal difference in enth- 
alpy potential. In this method, it is not possible 
to start with the coefficient and solve for the pre- 
dicted performance conditions except by trial and error. 

Hallett [23] with the help of CTI (Cooling Tower 
Institute) bulletins ATP-107B and ATP-127, presents 
methods which are used for calculating performance of 
a mechanical draft cooling tower. In this method, the 
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curves may be calculated using analytical methods 
similar to those used by the Cooling Tower Institute 
to determine performance levels from field test data. 

A format and calculation procedure for computing 

based on Tchebyshev's method for numerically evalua- 
ting the integral (eqn.2.5) as shown in detail in 
Chapter' III. This method is employed in the present 
work in evaluating the performance curves since this 
method gives consistent results over a wide variety of 
cooling ranges and wet bulb temperatures. The method 
also imparts itself to programming on a digital com- 
puter. 

CROSS F I- OW TOVRR PERFORMANCE : References on the cross 
flow principl e of water cooling originated in the 
early 1920’s. Some of the first cooling devices such 
as spray ponds and towers incorporated the cross flow 
principle in part, but the process was not analysed 
separately. 

The greatest advantages of the industrial cross 
flow tower are its design capability of water loadings 
to 50 m 3 /hr.m 2 of packing area (20 gpm/f t 2 of packing 
area) and its air velocities to '3.0 m/sec (600 f pm) . 

The methods used for estimating performance of 
counter flow towers cannot be applied to cross flow 
towers with the siame degree of accuracy, although 
there are certain similarities in the methods* The 
heat and mass transfer processes in both types of 
towers are based on the same potential for cooling. 
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The method of analysis and the process of integration 
are dependent upon the relative flow pattern of water 
and air in each type of tower. 

The cross flow tower involves a two-dimensional 
flow pattern in which water falls downward through the 
tower and the air is drawn horizontally through the 
packing. The enthalpy of the air changes not only in 
the vertical direction but also changes in the hori- 
zontal direction. 

In 1956, Zivi and Brand [24] presented a method 
for cross flow tower analysis based on the principle 
of enthalpy differential as the potential. They 
considered a vertical section through a cross flow 
cooling tower as shown in fig. (2,4). 


I* 


! 

1,1 

2,1 

3,1 

4,1 


t 

! 

Air Flow 
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2,2 

3,2 

4,2 
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j 

tA 

eg 

3,3 

4,3 

! 

J 



1,4 

2,4 

3,4 

4,4- 




Fig. (2.4) Vertical section through a cross flow 
cooling tower. 


The width of the tower is taken as unity. Hori- 
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zontal air flow is induced by the fan. The water 
enters at the top and flows uniformly downward through 
the fill. The positive X-direction is defined as the 
direction of the air flow, and the positive Z-direc- 
tion as the direction of the water flow. The water 
temperature entering the top of the fill is at the 
same temperature along the entire length. As the 
water flows downward, it is cooled. The water in the 
left hand portion of the fill is exposed to cooler and 
drier air than is the water in the right hand portion. 
Therefore, the rate of cooling of the water on the 
left is greater than that on the right. 

Thus, their method consists of a point- by-point 
determination of the water temperature and air enthal- 
py distribution by considering a small element of 
volume of the fill, and writing the equations of energy 
conservation and heat transfer. The contours of cons- 
tant average water temperatures are plotted on co- 
ordinates of non-dimensional tower height and depth. 

The dimensionless curves can be used to determine the 
characteristic required for a cross flow tower to 
meet a given cooling specification. 

Vouyoucalos [25] analysed the cross flow cooling 
tower by considering an element of volume with finite 
dimensions but sufficiently small so that the rela- 
tion h* 1 = f (T) could be approximated by a straight 
line. 

h' * = mT + p (2.12) 
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where , 


el = Constant, kcal/kg°C. 

p = Constant, kcal/kg. 

Considering enthalpy transfer equation; 

L C dT = K a (h" -h) dZ 
w 

2 


where. 


Z = 


L 

E 


f dT 
a ! h "-h 
l- 


Z = Side of the volume element, 
1,2 = Bottom and top of the tower. 

n 



: 




A.z 

~ i 

Ax 

i 

f 

i 



Pig. (2.5) lower co-ordinates (/^X = f\Z) 


(2.15) 

(2.14) 


in meters. 


Por this volume figure (2.5), equation (2.14) 
can he applied. By combining equation (2.14) with 
equation (2.12) and integrating 


(h” -h) 


Az K 


a 


(1 


,i , n~l 


l,n 

(h^TS)” 




- *3.n> I (2.15) 




- (h n -h) 




Since the variation of the group (h n -h) is not 
very great within a volume, the log-meaii can be approx- 
imated to the arithmetic mean 




48 


* 

« • 

2 (, - tj J 

(2.16) 




Considering the variation of h M within a 

volume 

element is not great 

- < h "> 3 ,n 


• 

« « 

K a 2(l i-l,n - T 1,n> 

( 2 . 17 ) 


L 11 "j, 11-1 “ h j,n-l - h j,n 

By 

(2.12), 



h 1 , = m T . , +p 

3,n-l 3-1, n v 

( 2 . 18 ) 

By 

basic equation 1 dT = G dh 



h 3 ,n = (T D-l,n " i + h 3 ,n-l 

( 2 . 19 ) 

By 

adding (2.17), ( 2 . 18 ) and (2.19) 



T . = B T . + C h . , + D 

3 ,.. 3 - 1 , n 3,n-l 

(2.20) 


where, 


3 = rows 
n = columns 


E a J\Z E a {\Z m 

r ^ ^ + i] 

B =s rr tc~<7 used in (2.20) 

E a & Z, 

[-2S— + l] 
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K a £z„ 

C t J 


C = 


[■ 


K _ a7S~Z” 


+ l] 


used in (2.20) 


pP % a ZL Z j 


D = 


sT a Az 


used in (2.20) 


2G- 


i3 


The enthalpy of the incoming air was h. 


the temperature of the incoming water was T 


o ,n * 


md 

both 


known. Equations (2j_9) and (2.20) can he used to 
find the enthalpy of the air ‘and temperature of the 
water in any part of the tower, starting from the 
volume element (0,0). 


This method could not be adopted in drawing the 
cross flow tower performance curves in the present 
work since, many variables were unknown and had to be 
obtained experimentally* 

Baker and Mart [19] applied the same analysis of 
Unit-Volume Coefficient' to cross flow tower although 
the calculations are more tedious. The mechanical 
integration was accomplished by dividing the cross 
section in a number of columns, each of which is sub- 
divided in a series of incremental volumes. The me- 
chanical integration of the coefficient for set of 
performance conditions is not too time consuming, but 
it is not possible to start with the coefficient and 
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solve for the predicted performance conditions except 
by trial and error. 

Halle tt [23] modified the Zivi and Brand method 
and showed the effect of variables which had been 
discussed but not completely defined by previous 
authors. The equations have been given in a general 
form. The performance curves are dimensionless . and 
are not related to a particular tower design. Hence, 
this method has been employed in the present work in 
drawing the performance curves of cross flow cooling 
tower. The format and calculation procedure to 
evaluate the enthalpy of air and temperature of water 
with the equations is shown in detail in chapter III. 
This method also contributes itself to programming on 
a digital computer. 



CHAPTER. III 


PERFORMANCE CURVES EVALUATION 

Various methods have been described in the pre- 
vious chapter to compute the performance curves for 
the mechanical draft cooling towers. In the present 
chapter, the methods proposed by Zivi and Brand [~4] 
and Halle tt [ 23 ] have been adopted for computing 
performance curves for the counter flow and cross flow 
towers, respectively. Computer programs have been 
developed using IBM 7044/1401 computers. The foll- 
owing industries have been taken into account as 
representative users of cooling towers: 

(l) Thermal power plants, 

(?) fertilizer plants and 
( 3 ) air conditioning plants. 

Tower performance has been studied for various 
design wet bulb and dry bulb temperatures for various 
major cities of India. The design values of wet bulb 
temperature and dry bulb temperature have been supp- 
lied by Paharpur Cooling Towers Pvt. Ltd., Calcutta as 
listed in table (3*1). 


U. L r - ; f'FUH 
CmWA L ! {BRAKY 

k 508SS 


An*. N5rn 
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Table { 3 . 1 ), Design wet bulb and dry bulb temperatures 
for various locations in India. 


Si. 

Ho. 

City 

-Lj ct « 

Deviation of 
I.S.T. from 
Solar noon, 
mins „ 

Design 

WBI, 

°C 

Des ign 
DBT, 

°C 

1 

Ahmedabad 

23 

+41 

’7.8 

41.1 

2 

Allahabad 

25 

+23 

28.3 

41.5 

3 

Amritsar 

32 

i-32 

28.3 

40.0 

4 

Asansol 

24 

-10 

28.3 

40.0 

5 

Bombay 

19 

+41 

27.8 

33.9 

6 

Calcutta 

23 

-22 

28.3 

36.7 

7 

Delhi 

29 

+23 

28.3 

40.5 

8 

Gauhati 

26 

-35 

28.3 

30.5 

9 

Gaya 

25 

- 8 

28.3 

4i. 7 

10 

Gwalior j 

26 

+ 19 

27.8 

41.6 

11 

Jamshedpur 

23 

-13 

28.3 

40.5 

12 

! J aipur 

27 

+29 

26.6 

40.0 

13 

Jodhpur 

26 

+39 

27.8 

41.5 

14 

Kanpur 

27 

+10 

28.3 

41.1 

15 

Lucknow 

27 

+ 8 

28.3 

40.0 

16 

Madras 

13 

+ 11 

28.3 

37.2 

17 

Nagpur 

21 

+15 

26.1 

42.6 

18 

Patna 

26 

- 9 

28.3 

39.4 

19 

Poona 

18 

+36 

24.5 

38.3 

20 

Trivendrum 

8 

+24 

27.2 

32.7 

21 

Yishakapa tnam 

18 

- 2 

28.9 

: 

36.6 


COUNlhR PLOW PERFORMANCE CUEVpS ; The curves for counter 
flow cooling tower are computed by using the Tchebyshev’s 
method proposed by the Cooling Tower Institute [26] ahd 
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by Hallett [23]. 


The calculation procedure is explained as follow- 
ing: 

In order to determine the cooling tower character- 
istic , it is necessary to integrate the definite 

integral (eqn. 2.5). 

> 

f dT 
I h”-h 

i 


or 



(3.1) 


between the known limits of water temperatures T-^ and 
T 2 . The Tchebyshev's method of integration is used 
according to which: 

/L 2 n 

,f(x) dx — 2Z f ( x ± ) (5.?) 

-1 J i=l 

where x^, indicate the re?l roots of the Tchebyshev 

quadrature polynomial for different values of n. For 
n = 4, the roots x^ are given as + 0. 187592^+ 0.2 and 

+ 0 . 794654$,+ 0.8. 
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Uow, in the integral (3.1), we can express 

T, + T 0 T, - To 
T = + ~ 6 


when 6=1, T = T ]L 
when 6 =-l, T = Tg 


T- 


T, 


Also dT = 


d6 


Substituting T in terms of 6 in the integral (3.1) and 
inserting the corresponding limits for 6, we gets 

, 1 

I. 


I = 


K a 


V -i 


1 T 1 


d6 


7J 



=( 




(3.3) 


Consider the height of a counter flow tower such 
that T 1 and T 2 represent the temperatures of the in- 
let and the outlet water, as shown in the figure (3.1). 
These temperatures correspond to the points + 1 and - 1, 
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respectively ox the interval (-1, + 1) considered in 
the int egral ( 3 .°) •' 



T r 0.l(T r T 2 ) 

T-L-0 .4(T 1 -T a ) 
(Il+T 2 )/2 

T 2 +O.4(0? 1 -T 2 ) 

t 2 +o.i(t -t 2 ) 


Fig. (3.1), nthalpy and temperature distribution in a 
counter flow tower. 


In order, now, to determine the temperatures at 
the four specified points in the interval (-l,+l), we 
use the relation 


T 


i 1+ 1 2 


when 6 =? -0.8, 

T = 0.5(T 1 +1 2 )-0.4(T 1 -T 2 ) 



) 6 


(3.4) 
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when 6 = -0,2, 

T = 0.5 (Tj+Tg) - 0.1 (T^-Tg) 

= T 2 + 0.4 (T^Tg) 

when 6 =+0.2, 

T = 1 1 - 0.4 (T^Tg) 

when 6 =+0.8, 

T = T x -0.1 (T 1 -T 2 ) 


(3.5) 

(3.6) 

(3.7) 


The relations (3.4 - 3.7) represent the tempera- 
tures at four different sections of a counter flow 
cooling tower. 


Similarly, we can express the enthalpy ass 


^ + h 2 (h 2 - h x ) 
h = ~ 2 ~ — ~ + 2 6 


(3.8) 


Following the same procedure, the enthalpy at the same 
corresponding sections of the tower can be obtained as: 


h = h x + 0,1 (hg-i^) (3.9) 

h = h x + 0.4 (hg-l^) (3.10) 

h = hg - 0,4 (hg-hi) (3.11) 

h = hg - 0.1 (hg-^) (3.12) 


These results are arranged in table (3.2) in a 
format used for our calculation procedure. 
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The design data on cooling towers for different 
services have been collected from various cooling 
tower installations in Kanpur city. They are shown 
in table ( 3 . 3 ), which consists of design values of 
water circulation rate (L), air flow rate (G-), inlet 
water temperature (T-^), outlet water temperature (T 0 ) , 
inlet wet bulb temperature (t^i) etc. 

Using the known data from table ( 3 . 3 ), and adopt- 
ing the procedure indicated in table (3*2'', the value 

17" IT 

of the tower characteristic — — is obtained. 

Example ; Consider a counter flow cooling tower inst- 
allation at Kanpur city for a thermal power plant of 
capacity 64 MW. The total heat to be dissipated is 

r 

128x10 kcal/hr. The inlet hot water temperature and 
outlet cold water temperatures are 43°C and 35°C, 
respectively (table 3 . 3 ). Hence, the cooling range 
is 8°C. The air flow rate is given as 8.66x10^ kg/ hr. 
The ambient air design wet bulb temperature for 
Kanpur is 28.3°C (table 3.1).' 

\ 

The following procedure is adopted to compute 
the design value of — . 

Q = L x A x IX T x C w 

where Q — Quantity of heat to be dissipated, kcal/hr. 
Al® = Cooling range, °C. 

A = Area of cross-section of the tower = 1100m 2 . 
C = Specific heat of water = 1 kcal/kg°C. 

“ h = ^TxAxC w 



Table (3.3), Design Data 
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s 
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water circulation rate L = 


128 x 10 6 ' 

(43-35) x 1100 x~T 

14,545 kg/hr. m 2 


Air flow rate G 


8.66 x 10 6 
1100 


Inlet water temperature T-^ 
Outlet water temperature 
.*. Cooling range (T^-Tg) 


7872 kg/hr. m 2 

1.8467 

43°C 

35°C 

8°C 

(i r i 2 ) | 

1.8467 x 8 

14.7735 k cal/ kg dry air. 


The enthalpy of air at the 
entering wet bull tempera- 
ture 28.3°C 


= 21.801 k cal/ kg dry air 

(from standard tables) . 
= "1B01 kcal/kg dry air 

= + (W e 

= 21.801 + 14.7735 


= 36*5745 kcal/kg dry air. 

USing the table (3*2), the design value of the tower 
characteristic, 5— = 0.9052* 


Corresponding to the design value of — | - for a 
specified data, there is a design point ? D for the 
parameters specified in table (3.3). This is shown 
in figure (3.2). For the Same value of the ^ ratio 
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and the range, there may be many combinations of the 
wet bulb and the cold water temperatures to give the 
same design value of the tower characteristic — ~ g , ~ » 
For each such combination, there should be a corres- 
ponding point in the chart. By joining all these 
points including the'hesign point*' we get, a curve 
which is termed as the Performance Curve . Similar 
performance curves are obtained for many other cool- 
ing ranges. The following procedure is adopted s 
Estimate a cold water temperature or approach for the 
wet bulb temperature being considered. Using table 
(3.2) , the value of is obtained for the estima- 


ted cold water temperature. 
K a V 


The calculated value of 
K a 7 


— is compared with the design value of — jj — . 

If it does not agree within a reasonable tolerance, 
use a trial and error solution and a refined estimate 
for the cold water temperature until the tolerance of 
— is met. The process is similarly repeated for 
different cold water temperatures. 


To avoid the tedious calculations by hand, a 
computer program has been developed (Appendix-A) for 
following computation ranges covering all possible 
cooling tower installations in India. 

(i) Wet hulb temperature : 15°C -- 30°C by an incre- 
ment of 1°C. 

(ii) Cold water temperature ; 15°C - 40°C by an incre- 
ment of 1°C. 

(iii) Cooling range : 1°C - 20°C by an increment of 1°C. 
The program reads the above input data together 
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■with, the constant 1/ G ratio and the design value of 

— jj — . A tolerance is given to get the accuracy of 
computation. The computation time required for cal- 

culating — ^ for the above ranges was 3-4 minutes 

on an IBM 7044/1401 computer. The subroutine ALPHA 
calculates the value of h M for a particular value of 
I, to be used in table (3.2). 

In plotting the curves, the wet bulb temperature 
is taken as abscissa and the cold water temperature 
as the ordinate. Cooling range was kept constant in 
drawing each curve, figures (3.2) - (3.9) show sets 
of performance curves for various cooling ranges and 
L/G ratios for different applications. 

Situation may, however, occur when it is nece- 
ssary to determine the tower characteristic for var- 
ious L/G ratios and approach values for a particular 
design wet bulb temperature and a fixed cooling range. 
The procedure adopted to draw the performance curves 
under these conditions is described below; 

In tbi • case, the cooling tower characteristic 
j L -jl J L 2 _g computed for a series of L/G ratios and cold 
water temperatures. In drawing a chart, the wet bulb 
temperature and the cooling range are kept constant. 

— is obtained by using the table (3.2). The foll- 
owing ranges were considered in order that the per- 
formance curves are applicable to the particular in- 
dustrial locations of the country. 



COLD WATER TEMPERATURE 


L/G =1*8467 
KaV/L =0'3053 
TOL =0*0 5 


30 ^ 


= Design Point 




•1 / 


/ y y 

y y y / 

/ ,x Y x 

x s s \/ / 


y V' 




/X 


X X / 

/ s 


y X X 


x / 



C APAC1T V= 64 MW 
L = 16.x IQ 6 Ka/hr 


HWT 

CWT 

WBT 

RANGE 


43 °C 
35°C 
28‘3 °C 
8°C 


' 20 25 28*3 30 

WET BULB TEMPERATURE °C 

FfG 3*2 COUNTER FLOW TOWER PERFORMANCE 
CURVES (THERMAL ROWER PLANT). 


WATER TEMPERATE! 


CAPACITY 

% 

L 

HWT 

CWT 

WBT 

RANGE 


= 64 MW 
= 16x1 0 6 
= 43 °C 
= 35 °C 
= 28'3°C 


Kg/hr 


L/G 

KqV 

L 

TOL 

P D 


= 2'6C 


010 

Des 

Poii 


20 25 

WET BULB TEMPERATURE °C 


28’3 
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WET BULB TEMPERATURE °C 
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(i) Wet bulb temperature; 24.5°C, 26.2°C, 27.8°C 
and 28.3°C. 

(ii) Cold -water temperature: 15°C - 40°C by an in- 
crement of 1°C , 

(iii) Cooling range: 1°C - 20°C by an increment of 1°C. 

( iv) L/G ratio: 1.0 — 3*0 by an increment of 0.4. 

A computer program has been written (Appendix-B) 
to evaluate the above values. The time required to 
compute these values was found to be the same as in 
the previous case. 

In plotting the curves, the L/C ratio is taken 
K a V 

as abscissa and — jj — as the ordinate. The wet bulb 
temperature and the cooling range are kept constant 
in drawing a particular set of performance curves. 
Figures (3.10) - (3.13) represent sets of performance 
curves for various cooling ranges and wet bulb temp- 
eratures . 

CROSS FL0¥__PER?0RMNCS CURVES : As already mentioned 
in the previous chapter, cross flow tower involves 
a two-dimensional flow pattern in which water falls 
downward through the tower and air is drawn horizon- 
tally through the packing. The enthalpy of the air 
changes not only in the vertical direction but also 
in the horizontal direction. 

Consider a vertical section through a cross flow 
cooling tower as shown in figure (2.4). The width of 
the tower is taken as unity. The positive X-direction 
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is defined as the direction of the air flew, and the 
positive Z -direction as the direction of water flow. 


Air 

Flow 



(1,D 

(2,1) 

(3,1) 

(4,1) 

7 

Az 

Jl_ 

j 

(1,2) 

(?,2) 

(3,2) 

(4,2) | 

i 


(1,3) 

(7,3) 

(3,3) 

( 4,3) 

(1,4) 

\ 

■ j 

(°,4) 

! (3,4) 

(4,4) 

j 


! 

i 



Fig. (2.4), Vertical section through a cross flow 
cooling tower. 


Assumptions : 

* 

(1) The hot water temperature entering the top of 
the tower is at the same temperature along the 
entire X~ boundary in the direction of air flow. 

(2) Entering wet bulb temperature of air is constant 
along the vertical Z-boundary. 

Consider a small element of volume of the fill, 
having dimensions /^X and j\ Z., respectively, . as 
shown in figure (3.14). The water temperatures enter- 
ing and leaving the element are T and T + A T , 
respectively. ,At I s found to be negative. Air 
enthalpies entering and leaving the element are H 

SI 
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and H a + / respectively, where A H a is positive. 

■ t 

•J/ water flow vL- 

r 


Air flow 



! 

IT 

k 

U-Ax_ i| 

I 


Pig, (3.14), Co-ordinates of a volume element. 


Therefore, energy balance for the incremental 
element is: 

<v 1 Ax [T„ - (T„ + Ay] = 0 AZ[(H a + 

AH a ) - H a ] (3.13) 

where, 

C = Unit heat capacity of water, kcal/kg °C. 

L = Water flow rate, kg/hr. m 2 horizontal area. 
T = Temperature of water, °C. 

G = Air flow rate, kg/hr.m 2 vertical area. 

H a = Enthalpy of air, kcal/kg dry air. 


or 


- = & Az As a 

G AH a 
E z ~ “ 1 EJ~ 


(3.14) 

(3.15) 
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If the distances /\X and J\.Z are reduced, approaching 
zero in the limit, eqn. (3.15) becomes 


~d’2 


w 


a 

i 


~d¥. 

Tk 


a 


(3.16) 


Equation (3.16) is the partial differential equation 
relating the rate of water temperature decrease with 
the rate of air enthalpy increase. 


When the transfer of heat takes place between 
the unsaturated air and a wetted surface (water sur- 
face) , the driving force is the "enthalpy potential". . 
It is defined as the difference between the enthalpy 
of air at the water temperature and that at the 
actual local wet bulb temperature. Therefore, using 
enthalpy potential as the driving force for heat 
transfer, the rate of heat transfer in the volume 
element is ; 

Ai = i c w Ax (-AV = k s Az 

(H w - V (3 - 17) 

where, 

/\ q = Heat transfer in the incremental volume 
A% A Z , k cal/ hr A 

K = Volumetric heat transfer coefficient 

cl 

for a specified cross flow packing, 
kcal/hr.m 3 (kcal/kg) . 


Volumetric heat transfer coefficient has been ex-* 
plained in Appendix-D. 
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= Enthalpy of saturated air at water temp- 
erature, kcal/kg dry air. 

The value of Z o is determined experimentally. It is a 
function of water velocity, air velocity, water temp- 
erature and .pa eking. 

Equation (3.17) can be written as; 

- 1 = K a < H v - H a ) < 3 - 18) 

If the distance /\ Z is reduced, approaching zero in 
the limit, equation (3.18) will be 


^w 

>7 



(H. 


w 


H ) 
a' 


(3.19) 


Equation (3.19) is the partial differential equation 
relating the rate of water temperature decrease to the 
magnitude of the local driving force, H^. - H . 

Combining equations (3.15) and (3.18), we have 


G y. Y 

Ah 


a 


a 


where , 


= K (H T - H ) 

a w a 

= < H w - Ha 5 

= hX (H w - H a ) 


(3.20) 


(3.21) 


Z^ /\X 

/\ X = -—q , represents a non-dimension- 


distance. 


If we consider two points * *1* ’ and ,, 2 t * separa- 
ted by the non-dimensional distance AX along the X- 
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boundary of the tower, the eqn. (3.21) may he express- 
ed as: 

A H a = zyt _5a*5>S] 

= [(H - H ). - (H - H ) ? 3 (5.22) 

2 LV w a 1 w a ^ 

Thus, the enthalpy difference for air is equal to the 
product of the average of the enthalpy potentials at 
the two successive horizontal points and the dimen- 
sionless distance "between them. 


In computing the air enthalpy for the points on 
the X — boundary, it is assumed that the hot water 
temperature entering the top of the tower is the same 
along the entire length in the direction of air flow. 


For a generalised solution of the air enthalpy 
along the X-boundary of the tower, consider the points 
along the X and Z boundaries as shown in figure (3.15) 

Applying the equation (3.22) for the points (i,j) 
and (i-1, j) in figure (3.15): 


H 


a(i, 3 ) ” H a(i--l,j) 


£X 


(H -H ) + (H -H ) | 

v w a w a 

(i-1, j) (i,3) . 

(3.23) 


A* 


H w(i~l, j) + B W(i,3) H a(i-1,3) H a(ir3) 


(3.24) 


a ( i , 3 * ' a{ i- 1 , 3 ' ) 

H y C i- j ) 


|_ Mi- 1 , 3 ) 
( i -1 s 0 ) 1 


Ax 


a(i-l, .j u H w( 1 - 1 , j ) +H w( i, j ) H a(i-l # j 


' 


i+1,3-1 ! i+5,D-l 


i+ 2,5 


i+3 , j 


i— l f j+1 i y< j+l i+l^j+1 i+2 y j+l i+3»3+l 

i-1 , j+2 _ i, j+2 i+1, j+2 i+2 , j+2 i+3, j+2 

i-1, a+3; i,j+3 i+1, j+3 i+2,j+3 i+3, 3+3 

- : 


Pig. (3.15), Two-dimensional arrangement of points in 


a cross flow tower. 


The equation (3.25) is used to determine the air 
enthalpy along the X-boundary. Since, it is assumed 
that the hot water temperature remains constant in the 
X-direction, the calculations are started with the 
follQwing values of i and 3 

i = 1 . 0 , 3 = 1 . 0 , /\± = 1 . 0 , Aj = 0.0 






* 
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Similarly, the water temperature is calculated for 
the points on Z - "boundary knowing that the entering 
wet bull temperature is constant along Z - boundary. 
Considering equation (3.18), 


(3.26) 

(3.27) * 

, - K a hZ 
where, / \ Z = — ^ - — 

In equation (3.27), A.Z is a non-dimensional tower 
length. Again, considering the points ' ’1* * and ’ ’2*’ 
in the vertical Z-direction, separated by a dimension- 
less distance /\,Z, the expression (3.27) may be writt- 
en as: 



(3.28) 

Thus, the enthalpy difference for water is equal 
to the product of the average of the enthalpy potent- 
ials at two successive vertical points, and the dimen . ! ion *, 
Idg-o distance A Z^ between those two successive points. 
Referring figure (3.15), and the points (i, 3-1) and 
( i,j) , the equation (3.28) becomes: 

■X In deriving this relation thb L.H.S. has been nulli- 
fied by the specific heat of water which is unity. 


Az 


(H 


w 


H a>l- < H W 


H ) 


a' 2 


At. 


w 


7yz~ 


E a {H w - V 

K a AZ 


1 < H w - V 


= - A'2 (H w - H„) 


a' 
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T w(i,j) T w(i,3-1) 


A 2 , v 

"2 (H w “ H a\. - -n 
^ ( i,3“D 


+ (H tt - H ) 

w a (i,3) 


(3.29) 


*** T w(i,j) T w(i,j-1) 


H w(i,o-l) + 


H w(i, 3 ‘) ” H a(i,j-l) H a(i,3)i 


(3.30) 


Since, it is assumed that the entering wet hulh 
temperature is constant in the Z direction, the cal- 
culations are started with the following values of 
i and 3* 

i = 1.0, 3 = 1.0, A 1 = °*°» A 3 = 1,0 

To solve equations (3.25) and (3.30), a computer 

program has been developed using IBM 7044/1401 

computers. The procedure to calculate and 

T /. is as follows; 
w(i, 3) 

X- -boundary 

The program begins by reading the initial hot water 
temperature entering the top of the- tower and the 
initial wet bulb temperature of air. These are de- 
noted by T„(i,;j-1) ond H a(1 _ lf3) . respectively. The 

program also reads the enthalpy of air at the enter- 
ing wet bulb temperature. ^a(i,3‘) aS So -^ ve< ^ ^ using 
the equation (3.25),, and having an increment = 1.0 

every time. 


i-3 In 


85 


-bound ary 

/ .x is calculated by using equation (3.30) as 
w(i,j) 

follows : 

(i) Assume (wet bulb temperature of the 

-.altering- air, °0 ) * 

(ii) Por the assumed value of T w (i 3 )* H w(i,j) iS 

calculated by using the standard tables and inter- 
polating results on the computer (Appendix-C^, 
subroutine ALPHA) . 

(iii) T w is computed by substituting lri 

equation (3.30). 

(iv) If 5? w ( i^ j) assumed ~ ^w( i, 3 ) calculated ’ 

tolerance, we stop and proceed for the next step 
to calculate -j) giving an increment 

A 3 = l.Q. A tolerance of 0.01 is assumed for 


accuracy. 

(v) If the tolerance is not met, 

_ T w(i, 3 ) assumed + T w(i, 3 ) calculated 

assume T / . .x = — 

w(.i,3J 2 

and the steps (ii), (iii), (iv) and (v) are repeated 
until the tolerance is met. 


In equations (3.25) and (3.30), a value of 0.06 
is taken for the dimensionless mesh size /\,X and /^Z. 

Int er ior P oint s 

By taking the known values of and from the 

preceding points, the values and ^ or 
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interior points, also can be obtained, again by using 
equations (5.25) and (5.30). Starting with i = 2 . 0 , 
j = 2 . 0 , both /\ i and /\ j are given an increment value 
of 1 . 0 , simultaneously, and the following steps are 
adopted: 

(i) Assume [H H / , 1 _ Tu TT 

^( 1 , 0 ) U a(i,3) J - [H »(i-l, 3 ) - H a(i-1,3) ] 

*• equation (3*25) becomes, 

H a(i,3) = H a(i“l, 3 ) + A* ( 3 . 31 ) 

ond, equation (3# 30) becomes, 

4 

AZ rrr 

- =1- GW-u-n - h. 


Iw (i,3) ^(1,3-1) 


2 ~ - "w(i,j-l) 
H w( i- 1 , 3 ) - H a (i-i, 3 )] 


a(i, j-l) 


(3.32) 


(ii) Compute and j) ^y using equations 

(3.31) and (3.32) 

(iii) For computed value of T , . H , . is cnl _ 

V(i,j) 1S caJ -~ 

culated by using the standard tables and inter- 
polating results on the computer (Appendix-^ 
and Appendix-^) . 

(iv) Calculate [H , . .* - h , ,7 

L w(i,j) a(i,j) •* 

(v) If [H , . .. _ h , ,7 ru 

assumed w(± f j)^ 

H 1 

a (i»j) J calculated, 

within the tolerance limit, we stop and proceed 

to the next step of calculation by giving incre- 
ments for A i - 1.0 and /\> j = 1.0. A tolerance 


of 0.01 is assumed for accuracy. 

(vi) If the tolerance is not met, take the average of 

the assumed and calculated values of fH , . - 

S n(± j) an(i t ^ e Bte P s (iii )t (i v) , (v) and 

(vi) are repeated until tolerance is met. 

For computing the values of H , H and T, a 

a w vt 

30 x 30 matrix of points has been considered with a 
mesh size of 0.06. The elements of the above matrix 
are divided into groups in a manner such that each 
element in a certain group gives the same value of the 
parameter to be determined. This classification is 
done to simplify the plotting of the performance cur- 
ves. The computer program developed is given in 
Appendix - C^. 

The values of unknown parameters are obtained 
similarly, for different inlet hot water temperatures 
and different entering wet bulb temperatures. 

These points are plotted to obtain a set of 
curves showing average hot water temperatures and 
enthalpies of moist air. These curves are dimension- 
less and are not related to a particular tower design* 
This permits the results to be applied for any cross 

flow toiler for which K , G- and 1 are known. 

a 

Figures (3.16) — (3.~?) show sets of performance 
curves for hot water temperatures of 43°C, 35.5°C and 
45.2°C and inlet design wet bulb temperatures of 24.5°C, 
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FIG. 3-21 CROSS FLOW COOLING TOWER 
PERFORMANCE CURVES. 
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26.2 C, 27.8 C and 28.3°C to suit various industrial 
applications for different locations of the country. 


CHAPTER IV 


RPSULIS AND Discussion 

die Tchcbyshev* s method for numerically , cvaluat 
ing -die I-Ic-rkol equations (2.5) end (2.5)' is consist- 
ent over a TT ido vanioty of cooling ranges mid wet 
bulb temperatures . As described in the previous 
cho.pt ir, this method has boon applied to determine 
the cooling to Tr or characteristic for various computa- 
tion ranges . Performance curves are dr run for those 
ranges both for counter flow and cross flow towers. 

Figures (3.2) - (3.9) show sets of curves for 
the counter flow tower, where the cold water tempera- 
ture has been plotted as a function of the design wet 
bulb temperature for different cooling ranges. The 
design conditions are shown by dotted linos. 

Four design wot bulb temperatures 24.5°C, 26.2°C 
27.8°C and 28.3°C have boon considered for various 
locations of the country Table (3.1). From tho per- 
formance curves, tho cold water temperature can be 
obtained at these wet bulb temperatures for a parti- 
cular cooling range depending on the nature of tho 
industry. The cooling range for various industries 
is as following [27] » 

3.5°C - 11°C for air conditioning and refri- 
geration systems 

5.5°G — 11°C for diesel engine cooling 

6.7°C - 11°C for steam surface condensers 


10 ? 


8.5 0 — 16 , 5°C for various industrial processes. 

The foiloxring examples Fill explain how the per- 
formance curves may he used by the manufacturer or the 
buyer. 


Cou nter _P1 ot-r Tower ; 

E xampl e 1 : It is desired to design a counter flow 
cooling tower for a thermal power plant in Kanpur 
where the design wet bulb temperature is 28.3°C. The 
following data are known for the thermal power plants 

Design dry bulb temperature at Kanpur = 41.5°C . 
Cooling range = 8°C . 

Height of the packing in the cooling tower, V ; 4 

meters . 

Quantity of heat to be dissipated = 100 x 10 6 kcal/hr. 

Consider a packing of rectangular slats as shown in 
figure (?.?). The rectangular slats are made of xrood. 
Choose water flow rate, L = 10,000 kg/hr. m 2 and air 
flow rate, G- = 55CO kg/hr. m 2 . Therefore, for the 
water to -air flow rate ratio = 1.8, the volumetric 
ma.ss transfer coefficient for the packing considered 
is given by, 

K^- = 2270 kg/m 3 „hr (kg of water/lcg of dry air), 
from table (2.1). 


* K a V _ 2270 x__4 

£ ~ lo", ooo 


0.901 


From the counter flow toxfer performance curves for the 
thermal poxrer plant figure (3.2), the outlet cold water 



temperature at a wet bulb temperature of 28.3 W C is 
= 35°C when ratio = 1.8 and — ^ = 0.901. Since, 

the cocling range is 8°C, the temperature of the in- - 
let hot ■’rat-v-'-r is 43 °C . The quantity of heat to be 
dissipated = 100 x 1C 5 kcal/hr. 

* • Area of the cooling tower cross section for a chosen 
water flow rate of 10,000 kg/hr, n* iss 


A = 


_0 ICO x 10 6 

x S~T x £ ~ 8.0 x~ 10, 000 

= 1250 m 2 


ilow, for the 
total number 


packing considered in figure 2.2), the 

of fill deck layers = 

w distance between two 

successive slats 


i. o 

= 08 = 

•• Total air flow rote required to maintain the out- 
let and inlet water temperature at 35°C and 43°C, 
respectively throughout the tower cross section, is- = 
5500 x 1250 
= 6.87 x 10 w kg/hr. 

Density Q f air at 41. 5° C DBT = 1.1255 kg/m 3 


Total volume air flow rate 


6.87 x 10 6 
1.1255 x 60 


m 3 / minute 


= 1.016 x 10 5 m 3 / minute 


Total H.P. 


required to drive the fans = 


1.016 x 10 5 „ 

_ — — S->s. 

226.53 6 


450.0 H.P. 


■£ A rule of thumb calculation for induced draft towers 
is that each 226.536 m 3 /minute of air exhausted 
requires one horsepower [3]. 


lo 1 ^- 

The cooling tower can be divided into many cells 
to ensure the continuous cold water supply even if one 
cell breaks down. Let us assume in this case that there 

4 

are four cells. Thus, four fans of about 120 H.P. each 
are required to be provided in the four cells so that 
the total H.P. = 450.00 as determined above. In order 
now to verify that the tower characteristic calculated 
in this example is the appropriate one, we use the per- 
formance curves given in figure ( 3 * 13 ), for the follow- 
ing data : W33T = 28*3% range = 8°C, approach = (35°C - 

28*3°C = 6. 7°C) , and ~ ratio = 1.8. 

G 

As it is clear from figures (3*2) - (3*9) j the wet 
bulb temperature has an effect on the outlet cold water 
temperature if the cooling range is kept constant* The 
cold water temperature increases with the increase of the 
wet bulb. It was found in the above example, that for a 
WBT = 28.3°C, the cold water temperature is = 3?°G. If 
the WBT* s are changed to 24*5°C, 26.2°C and 27*8°C which 
refer to different locations in the country, the outlet 
cold water temperatures will be 32.8°C, 33«8°C and 34.7°C 
respectively, for the same cooling range. 
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In figure (3*2) , the L/G ratio, r — and the pack- 
ing design remains the same, for different cooling ranges* 
Hence, for a constant tower size, the counter flow cooling 
toiler performs better at a wet bulb temperature of 28*3°C 
than at other wet bulb temperatures shorn above because the 
approach at this wet bulb is minimum = 6*7°c. Thus, it is 
the characteristic of a mechanical draft tower to give better 
performance at higher wet bulbs. 

Similarly, figures (3*3) - (3-9) and (3-10) - (3-12) 
may be used to determine the unkn own parameters and thereby, 
to design counter flow cooling towers for different applica- 
tions and at different locations. 

It is to be noted from figures (3- 10) - (3.13), that 

the approach to the wet bulb increases as the L/G ratio 

increases for a set of constant values of JLaJL, WBT and 

L ’ 

range . 

.Flow . Tower ; Figures (3*16) - (3*27) show performance 
curves for the cross flow cooling towers. The entering hot 
water temperatures of ^3°C, lf5*2°C and 35*5°C and the inlet 
air wet bulb temperatures of 24*5°C, 26«2°c, 27.8°C and 
28*3 C are considered. The curves indicate that for a parti- 
cular industrial application, the outlet cold water tempera- 
ture increases with the increase of the wet bulb temperature, 
if the temperature of the inlet hot water is kept constant. 

These curves are dimensionless and are not related to 
any particular tower design. This allows the results to 

be applied for any cross flow tower for which K a , G and L 
are known. 
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Tiia following formulae are used to calculate the 
dimensionless co-ordinates X and Z of various positions 
in the tower packing in the X raid Z directions, for 
known values of K o , G and 1. 




o '. 6 6 G 


( 3 . 21 ) 

(3.21) ’ 


where, 

X — Packing depth, in meters in the direction of 
air flow. 

The value 0.06 in the above equation appears as a re- 
sult of the mesh size chosen for the calculation. 
Similarly, 



(3.27) 


Z - 

Z = 0 ToTl (3.27)’ 


where, 

Z = Packing height in meters. 


In order to explain how the cross flow tower per- 
formance curves are used, we consider the sane example 
as stated for the counter flow tower case, whore the 
water flow rate, 1 = 10,000 kg/hr. m 2 and the air flow 
rate, G = 5500 kg/hr.m 2 . 


For the above water- to-air flow rate ratio - = 1.8, 
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the volumetric heat transfer coefficient, for the pack- 
ing being considered, is K r = 4500 kcal/hr.m 3 (kcal/kg 
, of dry nir) as explained in Appcndix-D. Let the pack- 
ing depth X, and the packing height Z bo 1.5 meters 
and 4.0 meters, respectively. 

Using relations (3.21) and (3.27)' , 

T - 4500 :: 1.5 
A ~ 0.0*6 z 5500 

= 20.0 

2_ 4_500 x 4.0 

= 0 X 6 'll 10000 

= 30.0 

Therefore, in figure (3.19), the proposed cooling tower 
would be represented by the rectangle bounded by the 
lines X = 0.0, X = 20.0 and Z = 0.0, Z = 30.0. The 
average cold water temperature shown at the po in t 
(X = 20.0, Z = 30.0) is about 33.3°C. This means that 
a tower hawing a dimensionless height Z of 30.0 and a 
dimensionless depth X of 20.0 would discharge water at 
about 33.3°C when operating at 28.3°C wet bulb temp- 
erature and 43°0 r>f inlet water temperature. Therefore 
the approach in this example is 5.0°G. The plan area 
of the tower is obtained as follows: 

Quantity of heat to be dissipated = 100 x 10 6 kcal/hr. 
Hot water temperature = 43°C. 

Cold water temperature= 33.3°C, from figure (3.19). 
Cooling range = (43.0°C - 33.3°C) 

= 9.7°C 

Water flow rate,! =10 ,000 kg/hr.m 2 
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Plan area of the cross _ q_ 2.00 x 10 ^ 

flow cooling tower c w ~ x ZTt x L 9.7 x 10,000 

= 1030 m 2 

Honce, comparing examples 1 and 2 we .find that for 
the same inlet hot water temperature of v3 C C, for the 
same amount of heat dissipation and for the same pack- 
ing height, the cold water temperatures for co un ter 
flow and cross flow cooling towers are 35°C and 33.3°C, 
respectively. Thus, a lower temperature is obtained 
in the cross flow cooling tower as compared to that in 
a counter flow installation. Also, the ground areas 
covered by the counter flow and cross flow cooling 
towers are 1250 m 2 and 1030 m 2 , respectively for simi- 
lar conditions. Hence, to get the same cold water 
temperature by the two types of cooling towers, the 
cross flaw tower requires a smaller area. There is a 
reduction of nearly 25 per cent in sj ze when compared 
to the counter flow cooling tower. This proves the 
well experienced, fact that considering all the fa,ctors, 
thvj cross flow "J overs are economical and give better 
performance than the counter flow towers , 

The cross flow tower performance curves a.re use- 
ful in predicting the unknown parameters for given 
design conditions and honce in designing the cross 
flow cooling towers. 


CONCLUSION : 

Performance curves have been drawn and analysed 
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both for the counter flow and cross flow cooling towers 
for the design conditions suitable for the thermo 1 
power plants , fertilizer plants and the air condition- 
ing plants at Kanpur. Only a particular cooling range 
of 8 C is chosen in drawing the performance curves for 

vs; _ in f lgures (3.10) - (3.13). 

More locations and cooling ranges have not been 
considered only to avoid the thesis becoming too volu- 
minous. The methods described herein and the computer 
programs developed thereby may, however, be used to 
obtain performance curves for any situation. An 
acceptable degree of accuracy in drawing these per- 
formance curves has been maintained. Tolerance of 
0.5 J C in approach and a tolerance of two to thrj > per 
cent in the value of — g--~ have been allowed. 

The numerical methods employed and the computer 
programs developed to obtain the performance curves 
have been explained. 

The present day cooling towers require careful 
designing and selection because of the increased costs 
and constraints on the availability of space. The 
performance curves developed in the present work should 
be of great use to the cooling tower manufacturers and 
buyers in the country in predicting the tower perform- 
ance and in its selection. 


llo 


REFERENCES 


1- H.C. A GRATAE , Evaporative Cooling and Cooling 
Equipment, Lecture Series I, Mech. Envv 
Dept., I.I.T. Kanpur, 1970. 

2. Charles A. Macaluso, Ecological Aspects of Cooling 

Systems, Cooling Towers, CEP Technical M^nu^l 
197-2, pp 111 - 114. ' ’ 

3. K.K. Mckelvey, Maxey Brooke, The Industrial Cooling 

Tower, Elsevier Publishing Company, L:nr!<~n, 


4. L.D. Berman, Evaporative Cooling of Circulating 
Water, Pergamon Press, London, 1961, p 182. 

i. illiam Sisson, Determining the * ize of Evaporation 
Ponds, Plant Engineering, Vol. 29, No . 1 . 5 ,' 

July 24, 1975, p 88. 

6. I.S. Murthy (Electrical Dept., Hyderabad, Indio), 

Natural Draught Cooling Towers, Ind. J. Power 
and River Valley Development, Vo. 21, No. 2, 
Pob. 1971, pp 57-60. 

7. G.B. Hill, Mechanical -Draft Cooling Towers, Chemical 

and Process Engineering, Heat Transfer Survey, 
Aug. 1959, pp 36-40. 

8. S.I. Malkin, Simplified Guide to Cooling Tower 

Economics, Plant Engineering, Vol. 25, No. 23, 
Nov. 11, 1971, p p 103-105. 

9. Ron Holzhauer, Industrial Cooling Towers, Plant 

engineering, Vol. 29, No. 15, July 24, 1975. 

10. D.R. Baker, H.A. Shryock, A Comprehensive Approach 

to the Analysis of Cooling Tower Performance, 
Trans. ASME, J. of Heat Transfer, Vol. 83, 1961, 
P 339. 

11. Lewis W.K., The Evaporation of a Liquid into a Gas, 

Trans. ASME, Vol. 44, 1922, p 325. 


Ill 


12 . 


13. 


14. 


15. 

16. 


17. 


18. 


19. 


20 . 

21 . 

22 . 


H... . Degler, Selection, Operation and Maintenance 
of Industrial Ceding equipment, Trans . ASHE. 
Vol. "'3, 1951, p 1031. 

Joseph Lichtenstein, Performance end Selection of 
Me chani cm 1-Dr of t Cooling lovers, Irons . AS HD 
Vol. 65, 1943, p 779. 

-4*-e. uondon, ■ .A. Mason, L.II.K. Boolter, Performance 
Chore ct eristics of a II echoni colly Induced Draft, 
Counter Plow, Pocked Coding Tower, Irons. AS ML . 
Vol. 62, 1940, pp 41-50. 

-T.J. Love, ond D.G-. Christie, Internotionol 

Development in Hoot Transfer, ASMB Port V. 
1961-62, pp 933-42. 

L.G-. Nnrny anakh edknr , IT. lTago.ro jo. Proceedings; 

Second national Heat and Mass Transfer, 
Conference, I.I.T. Kanpur, Dec. 13-15, 1973, 
p D— 27 . 

P.V. Bulanina, V.A. Morozov, A. Ye. Sukhov, High 
Capacity Counter Plow Cooling Towers, Heat 
Transfer, Sov. Res, Vol. 7, Ho. 3, May-June 
1975, pp 44-54. 

Robert Jorgensen, Pan Engineering, 6th ed., Buffalo 
Forge Company [1961], Buffalo, Hew York, 
p 226. 

D.R. Baker, L.T. Mart, Cooling Tower Characteristics 
as Determined by Unin Volume Coefficient, 
Refrigerating Engineering, Sept. 1952, pp 965- 


ASHRA? Guide and Data Book, Equipment, 1972, 
p 235. 

R.P. Moore, Hater Consumption Study for Hava jo 
Plant, AS CD J Power Div. , Vol.' 97, Ho.P04, 

Dec. 1971, Paper 8564, pp 749-66. 

Pfeiffer, ./hat Size of Cooling Towers, Chemical 
Engineering, Vol. 56, Ho. 4, 1949, p 98. 


112 


23. 


24. 

25. 

26. 
27. 


& * F * Ho.llett Pqv,^ 

c r os for 

% ss . for Power, Tel’. of 

" °* ’ cr * 1975, pp 503-- 

S . M . Z i vi "R ■r -n 

a. Xy x> 0 Jj m H rl A ■v-v nf 

Cooling lover, Refriv^r?+ S Cross 

-^g. 1956 , PP 31 ~ r ' tln 3 Engineering, ' 

S. Y.uyoucr.los, Cross Flou C'ni- 

wss'pplfe^ff’ VoL jSly yS0d ' 

° 00 ll ?fX°;' er Insti ^* Buuotm m-121, H0ust0n> 

S It GV G --’l^nlr^ / 7 i ri 

r-^nAncen Operation and 

lovers, PoFPr/ ™ 7 Re £° r V n Cooling 

P S.15. ' °^* 10 ^» JSTo. 3, March 1963 , 

Psy chrome try* tr-i^^?^' tion and 

’ ‘ n P-Xrai and Sons, Delhi. 


28. 


o o n o r> o o w o o r • 


w 


APPENDIX-* 


18 JOB 

IBFTC MAIN 


C 

c 


************** =M****« *************** ************** + ******* + ** + # + ^ 

THIS PROCECIRE IS TO DRAW THE PERFORMANCE CURVES SHOWING COLD 
WATER "EMPEPATURE AS A FUNCTICN OF WET EULE TEMPERATURE FOR VARIOUS 
RANG ES . 

ALG= WATE R TC AIR RATIO. 

HI =ENTHALPY OF ENTERING AIR iKCAL/KG DRY AIR. 

KAV/L=CQCL ING TOWER CHARECT P I ST I CS . 

RA=CCOL ING RANGE, CEG CENTIGRADE. 

TW=W£T BLLB TEMPERATURE OF THE ENTERING AIR, DEG CENTIGRADE. 
T2=TEMPEPATURE OF COLD WATER LEAVING, CEG CENTIGRADE. 

TOL = TOLERA NCE FOR KAV/L. 

VCON = DES IGN VALUE OF KAV/L. 


********************************** ****** *************************** ##:([))ei([i< 


DIMENSION T2 (50) *TW{ 50) * RA { 50 ),H1(50) 

COMMGN/CMD/R C 1 50 ) ,A(150J ,B( 150, NRR 

100 FORMAT (IX, 1311 1H* ) ) 

102 FORMAT!! X, 6C Cl H*) , * INPUT DATA *, 60 ( 1H* ) } 

106 FORMAT (8X» 10 1( 1H-J } 

200 FORMAT {8 FI 0.5) 

202 FORMATC10I3) 

204 FORMAT { /, 1>,*L/G =*,F8.4, 8X ,*DES IGN (KAVJ/L =*,F8.4,8X, 

l*TOL£RiNCE LIMIT =*,F8.4) 

206 FORMAT { /, 8X,10! (1H-) ,//, 

18X ,*wE~ BULB TEMP . * , 3X, *EN T ER ING AIR ENTHALPY*, 4X ,*COLD WATER TEMP 
2.*,13X,*RANGE*,18X,*KAV/l*,// ,8X,101 (1H-)) 

208 FORMAT ( 10 X, FlO .4 , 10X , F10 .4 »13X, F10 .4, 13X, F10.4,13X,F10.4) 

210 FORMAT I//, 1CX, *RANGE OF TEMPERATURE*, 15X,*A*,14X,*B*,/) 

212 FORMAT ( 5X,F10.4,4X,*TC*,4X,F10.4,5X,FI0.4,5X, F10.4,5X,F10.4I 

214 FORMAT (4F16.8) 

PRINT1 JO 
PRINT 10 2 

READ 2 >2 *N 1 , N2» NR »NRR 
PRINT2 j2^1,N2,NR,NRR 
PR INTI 00 

READ 200 ,( HI (K ),K=1,N1) 

PRINT200 ,l M(K),K = ;,N1) 

PR INT100 

READ 200 ,( TWIK ),K=1,N1) 

PRINT200 *(TW(K),K-1,M) 

PR INTI 00 

READ 200 »( 72 (K ) » K=1»N2I 
PRINT2 00,CT2CK),K=I,N2) 


-SSL 


o u u * 




3 

2 


PR INTI JO 

READ 2 JO,(RA(K>,K=’ -MR) 

PR INT2 j J ,lRfi (K ) » K= » MR ) 

PR IN T1 vO 
R { 1 )= 0 .0 
NN=NRR-i 
PRINT 21 j 
DO 4 K =1 ,NN 
L=K+1 

REAP21A,p( L),AlK),b(K) 

co R ^u E 12 ’ R,K) ’ R,L, ’ a,k, ' 8 ' k > 

PRINT 102 

READ 200 ,ALG,VCON, GL 
PR I N T2 ,i4 »A LG » V CON » T 0 L 
PRINT ; OC 

PRINT2J6 

DC 1 1=1 ,M1 
DO 2 J=i ,N2 
DO 3 K=i ,NR 
T1=T2( Jl+RAIK) 

H2 =H 1 { I ) +RA {K)*AtG 
TT=T2( J) 

RR =RA { K ) 

HH=H1(I) 

«ST=V D !^ T1,TT ’ hH ’ H2 - ALG ’ 8R * V ' 
IF (A8S (TEST ) .GT.TOL) GO TO 3 

?SwwuI’™ ,, ’ HU,, ' T2,J » »RACK) 

CONTINUE 

CONTINUE 

STOP 

END 


y , !4 


IBF7C DELTA 

SUBROUTINE DELTA^CALCULATE S^T(-** VA Tup *nc ***************** ************* 
1HE ' C «BVCHEFF FETHOOI KAV/LiTHE COOLING TONER CHARE 

**^* ***********************************************♦**♦**♦*****— #^—*^ 
suh=o“ TINE CELTA,T 'TT,HH,H2,ALG,RR,VI 


DC 1 K=1 ,4 

rrr=rr^alg 

GO TQ{ 2, 2,4,5) ,K 
** T =TT+0. i*RR 
HA =HH+ l*RRR 
GO TO :■ 

3 T =TT+;).4*gp 
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HA =HH+ > . 4* R R R 
GO TO 

4 T =Tl- <.4*RP 
HA=H 2 - J. 4*RRR 
GO TC 

5 T =T 1- 1*RR 
HA=H2- ;. 1*RRR 

6 CALL ALPHA (T »HW) 

HD = HW-H 1 

; SUR=SUH+:./HD 
V =RR*SLM *C .253 
RETURN 
ENC 

IBFTC ALPHA 

i* **$$*#****$*$#*** jJ*********************** $$#********## ****** *j| 

c THE SUBROUTINE ALPHA CALCULATES THE ENTHALPY HW FOR THE TEMP »TW * 

C A ANC B ARE CONSTANTS IN HW=A <TW) + B * 

C ***************************** ********************************************* 

SUBRGUT I NE ALPHAtT *HW 1 
CCMMCN/CMD/Rf 150) ,A(150) ,8115 Cl, NR 
DO 1 K=1 »NP 
IF(T.LE.RCK) } GO TO 2 
GO TC : 

2 J=K—1 

GO TC 4 
1 CONTINUE 
J=C 

4 HW=A(J ) * T+ E f J ) 

RETURN 

END 

ENTRY 


INPUT DATA 


16 36 24 91 


10 . G 4 

10.71 

11.42 

12.16 

12.92 

13.72 

14.55 

15.41 

16.32 

17.2 5 

18.23 

19.26 

20.33 

21 .45 

22.62 

23.83 

15.0 

16.0 

17.0 

18 . C 

19.0 

20.0 

21.0 

22.0 

23.0 

24.0 

25 .0 

26 . C 

27.0 

28.0 

29.0 

30.0 

15.0 

. 6.0 

17-0 

18 . C 

19.0 

20.0 

21.0 

22.0 

23.0 

24 .0 

25.0 

26 . C 

27.0 

28.0 

29.0 

30.0 

31.0 

32.0 

33.0 

34 . C 

35.0 

36.0 

37.0 

38.0 

39.0 

*• 0.0 

41.0 

42 . C 

43.0 

44.0 

45.0 

46.0 

47.0 

48.0 

49.0 

50.0 






0.0 i -0 

3.0 9.0 

16.0 .7.0 

2.0 

10.0 ; 

13.0 

R { K ) 

A ( K ) 

1.C000-J0CU 

0.^1000000 

2 ,cooouo.;o 

0.42399959 

3. COCOCOCO 

0.43559959 

4.C000 00 CO 

0.45159999 

5. COCOCOCQ 

0.^6100003 

6. COCOCOCO 

0 .-8000002 

7. COOOuOCO 

0.49159998 

8. COGQ'JOCD 

0.50699997 

9. 00 00 30 CO 

0.03100002 

lO.OOOOOOCO 

0.-^5000 02 

ll.OOCOuOCO 

0.56500000 

] 2. CO 00 00 CO 

0.58"C000^ 

13. coco 'o:o 

0.60959990 

i4.COOOCOCO 

0. -.2 9 00007 

15. COCOCOCO 

0.65699995 

16. COCOOOCO 

0.67000008 

17 . CQ0300C0 

0.70999992 

18. COCOOOCO 

0.74000001 

19. COCOCOCO 

0.75999999 

20. COOOCOCO 

o.eocoooo*’ 

21. CO GO ,t;o 

0. 825S9992 

22. 00 co ;cco 

0.86000001 

23. COCOCOCO 

0. 9055999” 1, 

24. COCO COCO 

0.93000007 

2 5. CO 00 JO CO 

0 • 98CCQ002 

26. COCO 30 CO 

l.u2999997 

27. COCOOOCO 

1.06999993 

28. CO 00 00 CO 

1 .12000012 

29. COCO 30 CO 

1.16599984 

30.C000D0C0 

1.21000004 

31. COCOOOCO 

1.27599997 

32. C00900C0 

1.34000015 

33. COOOuOCO 

1.39999986 

34. COCOOOCO 

1.47000003 

35. CO 00 00 CO 

1.52999997 

36.C000 ?0C0 

1.61000013 

37. COCOOOCO 

1.69000006 

38. COCOGO CO 

1.76999998 

39. COCOOOCO 

1. 86000013 

40. 00 00 00 CO 

1.93999958 


3.C 4.0 5.0 6.0 
11. G 12.0 13.0 14.0 
-9.0 2C.0 21.0 22.0 


B { K) 

2. 256CCCC1 
2.24200001 

2.21799999 
2.17000008 
2. 13399982 
2.03900003 
1. 967000G1 
1.86200047 
1.66959960 
1.543 99967 

1.34399986 

1.12359960 

0.8240013'. 

0.57659871 

0.18495947 

-0.01000023 

-0.64599962 

-1.16000175 

-1.51999855 

- 2.28300259 

-2.875S9725 
-3.50999832 
-4.61000061 
-5.C7000351 
-6. 270GC046 
-7.51999664 
-8.55999756 
-9.50559603 
-11.30999756 
-12.470C0122 

-14.56999969 

-16.43000031 

-18.34999847 

-20.66000366 

-22.69999695 

-25.50030000 

-28.38000488 

-31.33999634 

-34.76000977 

-37.87998962 


7.0 

15.0 

23.0 


41. CO 00 _ v OC0 
42.000000 CD 

43. COCOCOCO 

44. CO 00 JO CO 
45.00C0 iOCO 
46.0Q00C0C0 
47.0000 30 CO 

48. COCOCOCO 

49. COCOCO CO 
5 0.C000OG00 

51. COCOCOCO 

52. CO 00 00 UO 

53. CO COCO CO 

54. C000.j0C0 

55. COCOCOCO 
56. OOOOCO CO 

57 . COCOOOCO 

58. COCOOOCO 

59. COOOOOCO 

60. COCO 30 CO 

61. COOOOOCO 

62. COCOOOCO 

63. COOOOOCO 

64. COCOOOCO 

65. COCO CO CO 

66. COCO-COCO 

67. COCOCOCO 

68. COCO 33 CO 

69. COCOCOCO 

70 . CO 30 ; 3 CO 


2.35000019 

2.15000010 

2.15000000 

2.36999989 

2.49CC0025 

2.62999964 

2.76000023 

2.90999985 

3.05599994 

3.21999979 

3.41000080 
3.58999920 
3 .78000069 
4.CGOOOOOO 
4.22599954 
4.48000050 
4.73559977 
5.01959950 
5.34000015 
5.600000 38 

6 r» n r* f 

* ^ J W Joo 

6.39559962 

6.69559981 

7.30000114 

7.69559386 

8.30003114 

8.75999924 

9.35559962 

10.20000076 

10.79559924 


-42.27999878 
-46.38030498 
-50.58030183 
-55.73597498 
-6 1. 3 20 J 19 53 
-67.31997681 
-73.29598779 
-80.34597559 
-87.54998779 
-95.38558413 

-104. 89004517 
-114.06994629 
-123.95004272 
-135.60598535 
-148.02996826 
-161.78002930 
-176.33996582 
-192.29998779 
-210.86004639 
-226.20001221 

-250.2C3C1221 
-274.59997559 
-293.19995117 
-331. CCCCOOQO 
-356.59997559 
-395.60009766 
-428.55597559 
-468.79992676 
-523.20007324 
-564.59985352 


71. COCOOOCO 

72. COOOOOCO 

73. C000 30 CO 

74. COCOCO CO 

75. COOOOOCO 

76. COOOOOCO 

77. COCOOOCO 

78. COOOOOCO 

79. COOOOOCO 

80. COCOOOCO 


11.80000114 

12.59999847 

13.70000076 
14.89999962 
16.10000038 

17.59999847 
19.20000076 
21.10000229 
23.19999695 
25.60000229 


-634.60009766 
-691.39990234 
-770.60009766 
-858.19995117 
-947. C0024414 
-3059.49975586 
-1181.10009766 
-1327. 400'. 4648 
-1491.19970703 
-1680.80029297 


81. COCOOOCO 

82. COCOOOCO 

83. COOOOOCO 

84. COOOOOCO 

85. COOOOOCO 


28.59999847 

31.90000153 

35.70000076 
40.29999924 

45.70000076 


-1920.79980469 

-2188.09985352 

-2499.70019531 

-2881.49951172 

-2335.10009766 


ouuuuooooo 


ISA 
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IBJCB 

IB FTC MAIN 

********************************************************************** 
THIS PROGRAM IS TO COMPUTE THE VALUES FOR KAV/L AS A FUNCTION CF * 

L/C FOR VARIOUS APPROACHES ANC FOR A CONSTANT COOLING RANGE AND * 

A CONSTANT WET BULB TEMPERATURE. 

GL=VALLE OF L/G . 

HI =E NTH ALP Y OF ENTERING AIR f KCAL/KG DRY AIR. 

KAV/L=COCLING TOWER CHARECT P I ST ICS . 

RA-COOLI NG RANGE, DEG CENTIGRADE. 

TW=WET BLLB TEMPERATURE OF THE ENT ERING AIR, DEG CENTIGRADE. 

T2=TEMPE RATLRE OF COLD WATER LEAVING, DEG CENTIGRADE. 

3 ** 3= ************ =f * 4******************************************* 

DIMENSION T2(5D) ,TWC50) ,RA ISO ),H1( 50 ) ,GL(5Q) 

C0MM0N/CMD/R(150),AI 150) ,B ( 150 ) , NRR 
100 FORMAT ( 1 X, 131 ( 1H* ) ) 

102 FORMAT { 1 X, 60 (1H*) , *1 NPUT DATA*,60{ 1H* )) 

200 FO RM AT ( 3 FI C . 5 ) 

202 FORMAT (101 3) 

208 FORMAT { 2 X, 8 F 16 . 3 ) 

210 FORMAT (//,1CX,*RANGE OF TEMPEPATURE*,15X,*A*,14X,*B*,/> 

212 FORMAT ( 5X,F10.4 ,4X ,*T0* , AX ,F10 . 4 , 5X ,F 10 .4, 5X, F10.4, 5X , F 10 .4 ) 

214 FORMAT C4F16. 8) 

600 FORMAT (4JX,92(1H-) ) 

PRINTiOO 
PRINT 102 

READ 2 ;2 ,N1 ,N2 ,NR,NRR,NLG 
PR INT2~f2 »N 1 ,N2 »NR » NRR ,NLG 
PRINTIOO 

READ 200 ,(H1CK) ,K=1,N1) 

PRINT200,! H1(K),K=1,N1> 

PRINTIOO 

READ 2 00 ,( TW(K) ,K=1,N1) 

PR INT200 ,f T W CK ) ,K=1 , N1 ) 

PRINTIOO 

READ 200 ,( T2 (K ) ,K=1,N2) 

PRINT200 ,( T2(K ),K=1,N2) 

PR INTI 00 

READ 200 ,( RA(K) , K = 1 * NR) 

PRINT200,(RA(K),K=: ,NR) 

PR INT100 

READ 200f(GUK),K=l,NLG) 

PRINT2O0»(GL(K) ,K=1»NLG) 








R C 1 ) =0.3 
NN=NRR~: 

PR INTI 0 ') 

PRINT -.10 
DO 4 K=i »NN 
L s §<(•§* X 

READ214,R( L) ,A{K),SIK) 

PRINT ,;12» R(KJ ,R (L ) , AIK) ,B( K) 

A CONTINUE 
PRINT J Z'i 
DO 1 1 = 1 »M 1 
DO 2 J = ’.tNR 
DO 3 K =1 ,N2 
DO 5 L=i »NLG 
T1=T 2 1 K ) *R A I J ) 

H2 =H 1 1 I ) +R A ( J) #GL I L ) 

TT=T2( K ) 

RR=R A { J ) 

HH=Hll I ) 

GG=GL( L ) 

CALL DELTA IT1»TT »PF»H2»GG»RP,V) 

ASPRO=TT-TW II) 

PRINT 208, Hill ) »TW( I ) ,TT ,ASPRC,RA( J) ,GLIL) »V 

5 CONTINUE 
3 CONTINUE 
PRINT 600 
2 CONTINUE 

1 CONTINUE 
STOP 
END 

IBFTC DELTA 

C 4444 4444 44 4 4 44 4444 444444444! 44: 444 44 44444444 44444444444444^444444444 44444444 

C SUERCU-, IKE DELTA CALCULATES TEE VALUE CF KAV/L ,TFE COOLING TONER CHARECTE- 

C RISTICS BY THE '"CHEBYCHEFF RETHOC ) * 

l 44 4* 4444 4*44*4444444 44&44*# 4444444444444444444444444 444444 4444444444444444 

SUBROUTINE DELTA CT1* TT,HH,H2, ALG ,RR, V ) 

SUM=C. J 
DO 1 K =1 ,4 
PRP=RR*ALG 
GO TGC 2, 3 » 4 , 5 ) ,K 

2 T =TT+3.2*RR 
HA=HH+s . 1*RRR 
GO TC ;» 

3 T =TT+G.4*RR 
HA=HF+ J.4*RPR 
GO TO it 


O O O tJ 


A 


lit 


T =T 2- >.4*RR 
HA=H2- 1 * ^ RPR 
GQ TC 

5 T =T 1- • 1* R R 
HA =H 2- i . 1*RRR 

6 call alpha n»HW) 

HD =HW-H A 
, SU SUM+ 2. /HD 
V =RR*SlM*C.25) 

RETURN 

ENE 

I8FTC ALPHA 

THE SU8RCUT INE ALPHA CALCULATES THE ENTHALPY HW FOR THE TEMP »Tk * 

A ANC 8 ARE CONSTANTS IN HW=A (TW)+B * 

***************************** ********************************************* 


SUBROUTINE ALPHA ( T »HN) 

COYMCN/CMD/R (150) »A( 150 ) *B{ 150 »NR 
DO 1 K=1,NR 
IF(T.LE.RIK)) GO TC 2 
GO TC i 
2 J=K- 1 
GO TC 4 
1 CONTINUE 
J=0 

4 HW=A(JJ*T+B( JI 
RETURN 
ENC 

ENTRY ' 





INPUT 

DATA 




16 36 21 

91 16 







10.04 

0.71 

11.42 

12. 16 

12.92 

13.72 

14.55 

15.41 

15.32 

17.2 5 

18.23 

19.26 

20.33 

21.45 

22.62 

23.83 

15.0 

26 .0 

17.0 

18. € 

19.0 

20.0 

21. C 

22.0 

23.0 

24 .3 

25.0 

26.0 

27.0 

28.0 

29.0 

30.0 

15.0 

6.0 

17.0 

1E.C 

19.0 

20.0 

21-0 

22.0 

23. C 

24 .0 

25.0 

26. C 

27.0 

28.0 

29.0 

30.0 

31.0 

32 .0 

33.0 

34. C 

35.0 

36.0 

37.0 

38.0 

39.0 

40 .0 

41 .0 

42.0 

43.0 

44.0 

45.0 

46.0 


m 


4 .0 

-S .0 

49.0 

50 . C 

*» 





J.O 

1.3 

2.0 

2 .C 

4.0 

5.0 

6.0 

7.0 


?.o 

15.0 

9 .3 

: ?.o 

10.0 

18.0 

11. C 

: s.c 

12.0 

20.0 

13.0 

14.0 

15.0 


3.2 

0.4 

0.6 

€.8 

1.0 

1.2 

1.4 

1.6 


1.3 

2.0 

2.2 

2.4 

2.6 

2.8 

3.0 

3.2 


l FCP 

TEMPERATURE 

RANGES 

AND CCf'STANTS 

A AND 8 , 

SEE , 

APPENDI X-A 
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AFPENDIX-C 




IB JOB 

I BF "C MAIN 

**** 4 ***:<***** *************** ********************************************* 


u THIS PKOCEOLRE IS <Q CALCULATE THE DEMAND CURVES FOR VARIOUS CROSS-FLOW * 
C CONDITIONS. THE CROSS-FLOW COCLING TCwER INVOLVES A TWO-C IMENS IGNAl * 
C FL 06 PATTERN IN WHICH WATER FALLS DOWNWARD THROUGH THE TOWER AND AIR IS * 
C DRAWN HORIZONTALLY THROUGH THE PACKING. * 
C HA=ENThALP Y OF SATURATED AI R, KCAL/KG DRY AIR. * 
{ HE=E STIMATEC E. FHALPY, KCAL/KG . * 
L HX,HAX=M£SH SI ZE , D IMENS I ONU £S S , A VALUE CF 0.06 IS SELECTED. * 
C HW=ENTHALP Y OF SATURATED AIR AT WATER TEMPERATURE , KCAL/KG DRY AIR. * 
c HZ,HAZ=M;SH SI ZE, DIMENSIONLESS, a VALUE OF 0.06 IS SELECTED. * 
C TW=T EMP ERA TIRE OF WATER, DEG CENTIGRADE. * 
C TWB= WE T 2ULE T EMPE RATLRE OF THE ENTERING AIR. (DEG CENTIGRADE) * 
C TWE=ESTIMATED TEMPERATURE, DEG CENTIGRADE) * 


$ *** *********************** X* ************* ******************************** 


DIMENSION HA (30, 30) ,HW (30,3C ) ,TW( 30, 30 ) ,TW E (30,30) , HE( 30,30 ) , 

1HD(30,30 ), WOO, 30) 

COMMCN/CMD/R (150) , A ( 150 ) ,B ( 15C ) 

48 FORMAT ( 5 X, *LSES THE SECOND APPROXIMATION*) 

200 FORMAT (IX, 131 ( 1H* ) ) 

201 FORM AT ( 1 X, 60 (1H*),* INPUT DATA*,60( 1H*> ) 

210 FORMA! (//, 1CX,*RAMGE OF TEMPERATURE* , 15X ,* A*, 14X ,*B*,/ ) 

212 FORM AT ( 5X , FID .4 ,'->X ,*T0*, *X ,Fi 3.4, 5X,F10.4, 5X,FI0.4 ,5X ,F10.4 ) 

214 FORMAT ( 4 FI 6 . 8 ) 

215 FORMAT (3 F7. 3) 

PRINT2D0 

PR INT2 J 1 

NRR=SI 
ft( 1 ) =0. 3 
NN=NRR- ■ 

PRINT .1 : 

DO 4 K = l , N N 
L=K+1 

R£AD2i4,fi( L),A(K) ,B(K) 

PRINT Li 2, R ( K) ,R (L) ,A(K) ,B( K) 

^ CONTINUE 
PRINT2G0 

READ 2,5,HA(l,i) ,HW(1,1) ,TW(1,I) 

PRINT 215, HA(1,1),HW(1,1) ,T W ( 1 ,1 ) 

PRLNT2UU 
DO 30 1= ,3C 
30 HW(I ,1 ) = HW ( 1,1 ) 

DO 31 J=,,3G 




r> o o n 


31 HA ( 1 » J 1 (ifi. 1 
DO i 3 15=2,30 
,13 TW ( I E, . ) = 7 V ( i » 1 ) 

DO 1C 1 X =2 » 2 O 
I Y = I X- 
HX=0 .0 

10 HA (I X,. ) =HA (IY,1 ) + ( H X/2 . )*(HW (IY,1 ) + HK( IX, 1 )-HA ( IY, 1 > ) / ( 1 . +hX/2 . ) 

READ 2 5 *T W E 
PRIM! ^0 : 

D'J 1A IF=2 ,20 
IG= I F- 

TWE(1* IF )=(!*( 1, IG ) - ( TW { 1,1 G)-TWB) )/2. 

CALL ALPHA ( T WE ( 1 , IF) ,HW(1, IF) ) 

HZ =0.0-: 

14 TWU »IF) =TWC1, IG)-(HZ/2. )*( HW U ,IG )+HW( 1, I F)-HA(1 , IGl-HA ( 1, IF) ) 

18 DO 15 1H=2 ,20 

I F ( A BS (T W( 1,IHJ-7WE(1*IH) J .LE.0.01) GC TG 15 

PRINT 48 
GO TO ‘..6 

15 CONTINUE 
GO TC 19 

16 DO 17 IJ=2 ,3 0 
IA=I J~. 

TkEC 1 , 1 J )= ( TW ( 1 , I J )+ T WE( 1 , 1 J) )/2. 

CALL ALPHA (TWE (1 , I J) ,HW(1,I J) ) 

17 TW<1 ,1 J) =TW(1, IA)-(HZ/2. )*(FW(1»IA) + HW(1,I J )-HA ( 1 » I A )-HA( 1 » I J J I 

GO TC .8 

4 * *** 4 * 4 ******************* X* < 4 4 ** ** * *** * * * 4 = 4 = 4 :# 444444444444*************** 

TO CALCULATE THE VALUE AT THE INTERIOR POINTS KNOWING VALUES OF TW»HA,HW 
FROM PRECEDING POINTS 

4444 44444444444444444444444 ** $ 4 * 4 *** ********** **************************** 

19 DO 25 IT =2 ,30 

IU=I T- * 

DO 25 IV =2 ,30 
II=I V-l 

HE (IT, IV )=FW( IL, IV J-HA(IU,IV) 

FAX=Q. j6 
HAZ=C.Co 

HA (I T, IV 1= HA ( I U, IV )+(HAX/2» )*(HW( IU, IV ) -HA ( I U, I V )+HE ( IT , IV )) 

TW(IT, IV J=TW(IT,II )-(HAZ/2. )* (HW ( 1 1 , 1 1 )-HA ( IT » I I)+HE( IT, I V) > 

CALL ALPHA(TW( IT, IV) »8W( IT, IV } ) 

25 HD(IT,IV)=HW( IT, IV)-HA( IT, I V) 

40 DO 21 IP=2 ,30 
CO 21 10=2,3 0 

IF (AES (HC( IP, IQ)— HE( IP,IQ> ) .GT.0.01) GO TO 22 

21 CONTINUE 
GO TC 73 
DO 24 IK =2 ,30 


22 


# # 


O U U \J o o u u 


HEIIK.IL }= (HE( IK, ILJ+HOi IK, IL))/2. 

HA (IKt IL )= HA ( I M» I L ) + (HAX/2. )* (HVi( IM* ID-HA U M, I L) + HE ( IK, I L ) ) 

TW(I K, 1L )=TMIK, ID-tHAZ/2. )* IHW ( I K, I N )-HA ( IK , IN )+HE U K , I L > > 

CALL ALPHA (7W( IK, IL) »HW ( IK* IL ) ) 

2.A HC(IK, IL)=HW(IK»IL )-HA ( I K* I L ) 

GO TG * ■) 

23 CALL S 1 G HA C T rJ J 
PRINT2 > i 
CALL SIGNUHA) 

PR IN T2 u j 
6 CONT INUif 
STOP 
END 

IBFTC ALPHA 

***********:****************=***************=******************************** 
THE SUBROUTINE ALPHA CALCULATES THE ENTHALPY HW FOR THE TEMP»TW * 

A AND B ARE CONSTANTS IN HW=A (TW)+B * 

4 * 44 44 44 44 4444444 4444444444 44 4444444 444 *** ***** ************************ *** 

SUBROUTINE ALPHA (X,Y ) 

COMMON /C HD/RC150) » A ( 150) *B( 15Q) 

NR=9 1 

DO 1 K = 1 »MR 
IF(X.Lf.R( K))GO TG 2 
GO TC * 

Z J=K-1 

GO TO - 
1 CONTINUE 
J = C 

A Y=A{ J)*X+B { J) 

RETURN 
END 


IBFTC SIGMA 

44444444444444444444 *** 4**4 *4 4******************************************** 
THE SUBROUTINE SIGMA GROUPS THE REQUIRED ELEMENTS IN A 30*30 MATRIX 
ACCORDING TC TEMPERATURE ?W AND ENTHALPY HA) 

44 4*44:44: 444444444444444444 * 4 * 4444***************************************** 

SUBROUTINE S IGMA ( r W) 

DIMENSION TW (3'J»3C) 

TWMA =TW( i, 1 ) 

TWMI =TWC At 1 ) 

DO i K=1,3G 


# # 


DO 2 L = i *3 3 

I F IT W( K» L) .GT.TWMA ) TWMA=TW IK ,L ) 
IF(TWCK t L).LT.TWMI } TWMI =TW IK ,L) 

2 CGNT IN'IE 
l CONTINUE 

NI =1 FI X ( 7WMA ) 

N2 = IFI X ( TWMI ) 

N=N1 -N2+ i 
A = FLOAT(M Ml. 

DO 225 I 1=1. V 
A=A- 1. 

DO 2 27 KK=i ,30 
DO 2 26 L L= 1 * 30 
B=TW IKK, LL) 

EP SI =0 « i C 

IF (ABSIA-B) .GT.EPSI ) GO TO 22 6 
PRINT 22 fi» KK» LL » B 

226 CONTINUE 

227 CO NT INUE 
PRINT 40 u 

228 FORMAT (5 X, *TWI *,I2 12,* )=*,2X, F6.2) 
400 FQ RM AT 1 1 X, 3 C 1 1 H- ) ) 

225 CONTINUE 
RETURN 
EMC 

ENTRY 


INPUT DATA 

ENTER I EG AIR ENTHALPY, DEG CEN TRIGRADE 

ENTERING WATER TEMPERATURE , DEG CENTIGRADE 

WET BULB TEMPERATURE CF ENTERING AIR, DEG CENTIGRADE 


C FOR TEMPERATURE RANGES AND CONSTANTS A AND B, SEE APPEND IX-A * 
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APPENDIX - C 2 


Enthalpy of air rater vapour mixture at saturation [28] 

'0°C Dattm ' 


* 

jTemp. 

Enthalpy f h ) 

Temp . 

Enthalpy (h) 

T ' ' ' 

«■ Temp. 

iithalpy (h) 

(°c) 

(kcal/kg) 

( °C ) 

(kcal/kg) 

i (°C) 

i 

(kcal/kg) 

0 

° . 2 5 6 

30 

23.83 

60 

109.8 ! 

1 

2 * b 

1 31 

25.11 

61 

115.8 j 

2 

3.090 

1 32 

25.45 

62 

122.2 

3 

3.526 

! 33 

27.85 

63 

128.9 

4 

3.978 

34 

29.32 

64 

136.2 

5 

4.439 

35 

30.85 

65 

143 .9 

t 6 

4.919 

36 

32.46 

66 

152.2 

7 

5.411 

37 

34.15 

67 

161.0 

8 

5.918 

38 

35.92 

68 

170.4 

1 9 

6.449 

39 

37.78 

69 

180.6 

10 

6.994 

40 

39.72 

70 

191.4 

n 

7.559 

41 

41.77 

71 

203.2 

i 12 

8.144 

42 

43.92 

72 

215.8 

13 

8.754 

43 

46.17 

73 

: 229.5 

14 

1 9.383 

44 

48.54 

74 

244.4 

f 15 

10.04 

45 

51.03 

75 

260.5 

16 

10.71 

46 

53.66 

76 

278.1 

17 

11.42 

47 ! 

56.42 

77 

297.3 

18 

12.16 

48 1 

59.33 

! 78 

318.4 

19 

12.92 

49 

62.39 

79 

341.6 

20 

13.72 

50 1 

65.'61 

80 

367.2 

21 

14.5$ 

51 

69.02 

81 

395.8 

22 

15.41 

52 

72.61 

82 | 

427.7 

23 

16.32 

53 

76.39 

83 

463.4 

24 

17.25 

54 

80.39 

84 

503.7 

25 

18.25 

55 

84.62 

85 

549.4 

; 26 

19.26 

56 

89.10 

86 

602.3 

27 

20.33 

57 

93.84 

87 

663.4 

28 

21.45 

58 | 

98.86 

83 

734.7 

29 

22.62 

59 1 

104.2 

89 

819.1 


. . f 

V 

. . - . „ .j 

f 
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APPDHDIX— D 


Volumetric heat transfer coefficient [16]; ^hen it is 
very difficult to determine accurately the free sur- 
face of the liquid, e.g., on breaking up the flew of 
circulating writer into droplets, the volumetric of 
heat and mes-transf er coefficients are used, i.e., 
coefficients that are based, not on the unit surface 
cf the water, but on the unit active volume of the 

The volumetric heat transfer cceffic- 


cooling 

tewer. 
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Pig.(D-l), Comparison of performance characteristics 
of various packings (details of packing 
are given in table D-l) . 


lent is dependent upon so many factors like air 
flow rate, water flow rate, type of packing etc. 


me? 
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Various investigators have obtained the value for K 

y\ 

experimentally and they have been presented in figure 
(33-1). The details cf packings used by these invest- 
igators are given in table (I)~l). 


Table (D-l) 


SI. 

? 

! Investigator 

! -7-1 , 

i Packing 

| 

Horizontal 

jverticn 

Ho. 


1 

j . 

; S paeihg 

(Spacing 

1 

Lowe and Christie 

Corrugated 


f 



asbestos. 





cement louvres 

53 mm 

143 mm 

! 2 

Lichtenstein 

Vooden sl.ats 



♦ 


(10x8) mn 

j 10 mm 

380 mm 

3 

U chi da et al 

Honey comb j 





board, plus t it. - j 




! 

i 

j 

i j 

[ zaticn paper j 

112 mm 

65 mm 

4 

Simpson et al 

i 1 

| Masonite sheet j 

15 mm 

- ■ 

5 

Uchida and 

Inclined wooder 

t 


i 

Shnrah 

slats (10 mm) j 

72 mm 

67 mm 
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